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identical  conditions  hav^  a  lattice  constant"" af  5.6532A,  which 
matches  GaAs  to  within  (3  C  lO^^A^^flck  reflection  Laue  patterns 
show  that  the  {110}  layers  are  chalcopyrite ,  while  the  {111} As 
layers  are  sphalerite.  The  Laue  patterns  also  indicate  that  the 
{110}  epitaxy  is  mixed  with  the  c-axis  out  of  the  growth  plane v 
A  possible  explanation  for  these  results  is  that  the  {110}G« 
surfaces  in  liquid  phase  epitaxy  reconstruct  into  chalcgpyflte-like 
chains  in  the  [001]  direction,  while  the  {111} As  su££d6es  do  not. 

This  would  promote  the  formation  of  chalcopyri^e--^nSnP2  with  anti¬ 
phase  boundaries  on  {110}  surfaces. 

Using  back  reflection  Laue  measuawttents  we  have  examined  the 
substrate  orientation  i]j,m  mli  in  i  nl  antiphase  domain  boundaries 

on  GaAs  by  liquid  phase  epitaxy.  These 
results  together  with  those  of  other  workers  indicate  that  {211} 
substrates  are  best  for  the  growth  of  chalcopyrite  layers  with  c=2a, 
while  {100}  substrates  are  best  for  those  chalcopyrites  with  c<2a. 

The  suppression  of  antiphase  domain  boundaries  can  be  related  to 
improved  electrical  properties  for  the  chalcopyrite  layers.  In 
general,  crystal  orientations  with  large  c-axis  components  parallel 
to  the  growth  plane  tend  to  be  suppressed.  This  suggests  that 
differences  in  lateral  growth  rate  from  nuclei  with  different 
orientations,  rather  than  site  preference,  is  a  dominant  factor  in 
the  suppression  of  antiphase  domain  boundaries. 

Misordering  defects  which  occur  in  the  growth  of  chalcopyrite/ 
sphalerite  systems  have  been  investigated  through  the  LPE  growth 
of  ZnSnP2  on  GaAs.  X-ray  diffraction  measurements  reveal  the 
presence  of  misordering  defects  associated  with  multiple  tetragonal 
orientations  in  the  epitaxy.  Growth  on  {211}  oriented  substrates 
reveals  the  lowest  concentrations  of  orientation-related  defects 
while  growth  on  vicinal  {110}  oriented  substrates  indicates  the 
strongest  cation  site  differentiation  required  to  produce  long-range 
ordering.  The  use  of  critically  misoriented  {110}  substrates  is 
proposed  for  the  growth  of  highly  ordered  chalcopyrite  layers  on 
sphalerite  substrates. 

The  mixing  enthalpy  of  ternary  tetrahedral  semiconductor 
alloys  is  fairly  well  described  by  regular  solution  theory,  with 
a  thermodynamic  interaction  parameter  that  is  sensitive  to  the 
lattice  spacing  of  the  binary  constituents.  We  derive  an  estimate 
of  the  interaction  parameter  from  a  model  which  ascribes  the  mixing 
enthalpy  to  bond  distortions  associated  with  the  alloy  formation, 
and  relates  these  to  the  macroscopic  elastic  properties  of  the  crystal. 
Numerical  estimates  are  given  for  the  18  alloys  with  cations  Al, 

Ga,  In  and  anions  P,  As,  Sb  and  are  compared  with  experimental 
values  and  alternative  models.  To  within  a  single  adjustable 
parameter,  the  predictions  agree  with  experiment  and  are  consistent 
with  those  of  the  delta  lattice  parameter  (DLP )  model.  A  further 
calculation  of  the  elastic  energy  associated  with  composition 
fluctuations  (clustering)  in  these  alloys  indicates  that  this  energy 
is  sufficient  to  suppress  clustering  above  the  critical  mixing 
temperature. 
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20.  Abstract  (Continued) 

Spatial  composition  fluctuations  in  semiconductor  alloys, 

AxBi_xC  involve  elastic  energy  when  the  lattice  parameters  of 
AC  and  BC  differ.  Therefore,  the  fluctuations  are  expected  to 
be  reduced  below  those  exhibited  by  a  purely  random  distribution 
of  the  constituent  ions.  We  estimate  the  extent  of  the  reduction, 
and  relate  it  to  recent  observations  of  Raman  linewidths  and  line 
shapes  in  GaxIni_xAs  and  GaxAli-xAs. 

Below  band-gap  photovoltaic  response  with  a  half-power  point 
at  1.31eV  is  observed  in  p-type  Ini-xGaxP  on  n-type  GaAs  hetero¬ 
junctions  with  {111}  interfaces.  This  response  is  apparently  due 
to  photon-assisted  tunneling  of  carriers  across  the  interfacial 
energy  gap,  which  is  smaller  than  the  energy  gaps  of  the  constituent 
materials.  The  same  mechanism  can  be  employed  in  other  hetero¬ 
junctions  with  appropriate  energy-band  lineups  to  obtain  longer 
wavelength  infrared  response. 

Below  bandgap  emission  and  absorption  are  observed  in  p-ZnSnP2 
on  n-GaAs  hetero junctions.  The  band  lineup  of  this  heterostructure 
system  suggests  that  the  below  bandgap  transition  is  between  the 
valence  bands  of  GaAs  and  the  conduction  band  of  ZnSnP2  at  the 
interface.  The  predicted  value  for  the  interfacial  gap  of  Ej=1.27eV 
compares  with  an  experimental  value  of  1.3 lev 

Advances  in  technology  have  made  possible  the  fabrication  of 
rapidly  varying  heterostructures  which  hold  the  promise  of  important 
applications.  We  develop  a  set  of  approximate  treatments  of  electron 
states  in  a  variety  of  layered  heterostructures.  The  approximations 
are  all  based  on  the  concept  of  one-band  generalized  Wannier  functions. 
Following  a  discussion  of  the  validity  of  this  representation,  we 
apply  it  to  an  evaluation  of  the  bound  states  in  a  narrow  quantum 
well  in  GaAs ,  which  clearly  demonstrates  the  mixing  of  main  and 
satellite  valley  states  as  well  as  the  contribution  of  evanescent 
states,  and  of  the  states  of  a  super lattice  in  a  model  structure 
of  up  to  20  quantum  wells.  As  a  final  example  we  discuss  the  appli¬ 
cation  of  generalized  Wannier  functions  to  the  matching  of  electronic 
states  at  a  hetero junction  between  two  model  band  structures  with 
different  effective  masses,  and  compare  the  formalism  with  alterna¬ 
tive  approaches  to  this  problem. 

Recent  III-V  alloy  formation  models  indicate  that  negative 
charge  transfer  energy  can  overcome  positive  bond  distortion 
energy  to  stabilize  long-range  order  or  compound  formation. 

Although  a  report  of  a  layered  ordering  in  AlxGai_xAs  tends  to 
confirm  this  result,  we  have  as  yet  obtained  no  convincing  evidence 
for  such  ordering  in  InxGa]__xP.  Compound  formation  in  these  alloys 
may  be  critically  dependent  on  growth  temperature  and  kinetic  factors. 
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1.  INTRODUCTION 

Although  there  is  substantial  interest  in  III-V  semiconductor 
alloys  for  electronic  and  optoelectronic  devices,  relatively 
little  is  known  about  the  distribution  of  atoms  in  these  materials 
or  the  effects  of  non-random  distributions  on  device  performance. 

In  many  of  these  alloys,  however,  non-random  distributions  of  atoms 
are  expected  to  be  important  from  thermodynamic  considerations. 

This  is  particularly  true  at  the  low  temperatures  commonly  used 
for  epitaxial  growth.  The  available  evidence  suggests  that  short- 
range  clustering  of  like  atoms  or  short-range  ordering  of  unlike 
atoms  can  produce  device  problems  such  as  excess  noise  and  leakage 
current,  premature  voltage  breakdown,  and  lower  carrier  mobilities. 
Long-range  ordering  of  unlike  atoms,  however,  could  potentially 
yield  III-V  ternary  compounds  with  properties  superior  to  their 
parent  alloys.  In  epitaxial  growth,  substrate  effects  such  as 
lattice  match  and  non-equivalent  sublattice  sites  are  expected 
to  have  a  strong  influence  on  these  phenomena.  Also,  two  of  the 
important  III-V  alloys,  Gaxln^  P  and  GaxIni-xAs»  can  be  epitaxially 
grown  matched  to  GaAs  and  InP,  respectively,  at  compositions  near 
the  point  of  greatest  ordering  probability.  The  objective  of 
this  work  is  to  investigate  various  aspects  of  clustering  and 
ordering  in  III-V  alloys,  including  the  interfacial  properties 
of  the  materials. 
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Liquid  Phase  Epitaxial  Growth  of  ZnSnP2  on  GaAs 

G.  A.  Davis  and  C.  M.  Wolfe* 

Semiconductor  Research  Laboratory,  Washington  University,  St.  Louis,  Missouri  6 3130 


ABSTRACT 

ZnSnPi  is  a  potentially  useful  semiconductor  which  can  have  either  the  sphalerite  structure  or  the  chalcopyrite  struc¬ 
ture  with  no  tetragonal  distortion  (c  -  2a).  We  have  grown  ZnSnPt  on  various  orientations  of  GaAs  by  liquid  phase  epitaxy 
and  found  that  the  best  growth  occurs  on  (1 10)  surfaces.  Double-crystal  x-ray  diffraction  measurements  indicate  that  these 
{110}  layers  have  a  lattice  constant  of  S.6S07A,  while  {1 1 1}  As  layers  grown  under  identical  conditions  have  a  lattice  constant 
of  5.6S32A,  which  matches  GaAs  to  within  ±2  x  10'4A.  Back  reflection  Laue  patterns  show  that  the  {110}  layers  are 
chalcopyrite,  while  the  {1 11}  As  layers  are  sphalerite.  The  Laue  patterns  also  indicate  that  the  {110}  epitaxy  is  mixed  with  the 
c-axis  out  of  the  growth  plane.  A  possible  explanation  for  these  results  is  that  the  {1 10}  GaAs  surfaces  in  liquid  phase  epitaxy 
reconstruct  into  chalcopyrite-like  chains  in  the  [001]  direction,  while  the  {111}  As  surfaces  do  not  This  would  promote  the 
formation  of  chalcopyrite  ZnSnP,  with  antiphase  boundaries  on  {110}  surfaces. 


ZnSnPt  is  an  interesting  semiconductor  for  the  study 
of  chalcopyrite  growth  on  sphalerite  GaAs.  The  prim¬ 
ary  reason  for  studying  this  system  is  the  similarity 
between  the  structures  of  GaAs  and  ZnSnPi  (1).  The 
a  to  a  mismatch  of  ZnSnPj  on  GaAs  is  small  which 
allows  for  single  crystal  liquid  phase  growth.  In  addi¬ 
tion,  ZnSnPt  exhibits  no  tetragonal  distortion  so  the 
c  to  2a  mismatch  is  the  same  as  the  a  to  a  mismatch. 
Thus,  ZnSnPt  grown  on  any  substrate  orientation  will 
not  incorporate  any  c  to  2a  mismatch.  Since  most 
chakopyrites  exhibit  some  tetragonal  distortion,  epi¬ 
taxial  growth  is  generally  restricted  to  {100}  oriented 
substrates  to  avoid  c  to  2a  mismatch. 

ZnSnPt  is  also  of  interest  because  solution-grown 
platelets  exhibit  both  chalcopyrite  and  disordered 
sphalerite  phases  apparently  having  the  same  lattice 
constant  (1,  2).  The  sphalerite  phase  dominates  in 
platelets  grown  using  large  cooling  rates  (~  50" C/hr), 
while  the  chalcopyrite  phase  dominates  for  smaller 
cooling  rates  (~  5*  C/hr) .  The  existence  of  these  two 
phases  allows  us  to  evaluate  the  growth  conditions 
leading  to  chalcopyrite  growth  and  provides  insight 
for  other  chalcopyrite-on-sphalerite  systems. 

Growth  Procedure 

Liquid  phase  epitaxial  growth  of  ZnSnPs  on  GaAs 
was  performed  in  the  open  tube  system  shown  sche¬ 
matically  in  Fig.  1.  The  starting  materials  used  in 
this  system  were  90.999+%  pure  Zn  and  Sn  and  SnF* 
which  was  synthesized  from  99.999+%  pure  Sn  and 
red  P«  in  a  sealed  tube  process  (3).  The  graphite  boat 
shown  in  Fig.  1  was  designed  for  growth  of  materials 
for  which  there  is  little  phase  diagram  information. 
During  homogenization  of  the  growth  melt,  excess 
SnPj  floats  to  the  top  and  maintains  a  saturated  solu¬ 
tion  as  phosphorus  evaporates  into  the  gas  stream. 
During  growth,  the  homogeneous  solution  contained 
in  the  slider  is  positioned  over  the  substrate  and  the 
temperature  of  the  furnace  is  reduced  linearly.  Since 
little  phosphorus  loss  to  the  gas  stream  was  observed 
for  the  growth  conditions  used,  the  composition  of 
the  melt  was  adjusted  so  that  there  was  little  or  no 
excess  SnPj  in  the  melt  after  saturation.  The  equiva¬ 
lent  melt  compositions  used  for  various  growth  tem¬ 
peratures  are  shown  in  the  partial  phase  diagram  of 
Fig.  2.  In  all  growth  runs,  the  amounts  of  Zn,  Sn,  and 
SnP*  were  adjusted  to  give  an  equivalent  ZnSnPj 
concentration  in  Sn  with  no  excess  Zn  or  P. 

The  cooling  rates  used  in  these  growth  experiments 
ranged  from  1*  to  5"  C/hr.  Cooling  rates  greater  than 
3 •  C/hr  were  not  investigated  due  to  the  tendency  of 
ZnSnPi  to  grow  in  the  sphalerite  phase  at  these  rates. 
A  homogenization  period  of  2  hr  and  growth  periods 
of  8-18  hr  were  used.  The  homogenization  or  initial 
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growth  temperature  ranged  from  530*  to  840*  C.  A 
hydrogen  flow  rate  of  30  ml/ min  was  maintained  in 
the  furnace  tube  during  homogenization  and  growth 
for  all  experiments. 

The  GaAs  substrates  were  carefully  prepared  by 
lapping  and  chemical-mechanical  polishing  both  the 
front  and  back  surfaces.  After  scribing  and  breaking 
the  substrates  to  the  proper  dimensions  for  the  boat 
(0.25  x  0.25  in.),  they  were  cleaned  in  boiling  tri¬ 
chloroethylene,  acetone,  and  methyl  alcohol.  A  5:1:1 
or  10:1:1  H*SOe:  HyO*:  HyO  chemical  etch  was  used 
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R».  I.  Schematic  illestrohon  of  the  lit  host  used  to  f row 
ZnSaPy  am  GaAt. 


Fip.  1  Initial  (rawth  tempera  tare  at  a  f— cries  et  ZnSnPi 
concentre  ties  far  sotaratad  er  sUphtty  sops  tee  toretsd  satetians. 


Vol  130,  No.  8 


EPITAXIAL  GROWTH  OF  ZnSnP, 


1409 


prior  to  loading  each  substrate  into  the  boat.  This 
etch  removed  at  least  10  pm  of  material  from  the 
growth  surface.  The  10:1:1  etch  used  in  later  ex¬ 
periments  caused  less  rounding  of  the  surface  and  pro¬ 
duced  more  uniform  growth. 

Layer  Morphology 

Initial  growth  experiments  were  performed  on  {100}, 
{110},  {111},  and  {211}  oriented  Cr-doped  substrates 
to  determine  the  optimum  orientation  for  growth. 
Typical  layers  grown  on  some  of  these  orientations 
are  shown  in  Fig.  3, 4,  and  5. 

{100},  {211}  Ga,  {211}  As,  and  {111}  Ga  oriented 
substrates  produced  growth  which  was  quite  rough 
and  unsuitable  for  electrical  measurements.  Growth 
on  {111}  As  oriented  substrates  had  rough  interfaces 
but  uniform  surfaces  with  very  good  morphologies. 
This  is  consistent  with  observations  of  predominant 
{111}  facets  in  solution  grown  platelets  (1)  and  with 
experiments  in  which  we  have  observed  low  growth 
rates  in  the  <H1>  directions.  Layers  grown  on  {110} 
oriented  substrates  revealed  much  smoother  interfaces 
but  had  surfaces  with  more  structure  and  tended  to  be 
of  nonuniform  thickness.  Numerous  growth  experi¬ 
ments  were  performed  on  {111}  As  oriented  substrates 
to  reduce  the  amount  of  interfacial  roughness,  but 
little  improvement  was  noted.  In  later  experiments 
performed  using  {110}  oriented  substrates,  we  con¬ 
centrated  on  improving  the  uniformity  of  the  layers 
and  eliminating  large  growth  structures.  By  chang¬ 
ing  the  final  etch  and  adjusting  the  growth  param¬ 
eters,  substantial  improvements  in  the  {110}  oriented 
growth  were  noted. 

For  {111}  As  oriented  substrates,  the  best  growth 
eras  obtained  using  an  initial  growth  temperature  of 
M0*C  with  a  cooling  rate  of  l*C/hr  The  melts  for 
these  growth  runs  had  an  equivalent  ZnSnPi  concen¬ 
tration  of  1.4-1.9%  with  one  gram  of  Sn  used  per 
run.  For  initial  growth  temperatures  below  530* C,  no 
layers  could  be  grown,  while  for  initial  temperatures 
above  600*C.  the  morphologies  and  electrical  charac¬ 
teristics  of  the  layers  degraded.  The  grown  layers 
were  generally  3-19  pm  thick  except  for  those  grown 
at  high  temperatures,  which  were  up  to  30  pm  thick. 
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PI*.  X  (a.  Me)  Surface  aag  (b,  bsttm)  crow  wctieaal  vim  of 
typical  ZaSePt  ptowtb  oa  {100}  srioaMg  GaAs. 
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Of.  4.  (s,  toe)  Sarfaco  aag  (b,  bettors)  crow  sactMael  view*  of 
typical  ZaSaPp  leyots  |im>  oa  {111 }  As  arise  tog  GaAs.  Tbs  stoic - 
taro  of  tbs  sarfaco  vise  was  sahoacag  strep  s  Noras raki  istsr- 
fa  roars  tor. 


100pm 


5pm 


B*.  5.  (a.  Me)  Sarfaca  mi  (b,  bettors)  crow  tscfiesol  vim  of 
a  typical  ZaSaPt  leyor  prove  oa  {110}  arise  tori  GaAs. 

The  conditions  producing  good  growth  on  {110} 
oriented  substrates  are  more  restrictive  than  those 
for  {111}  As  growth.  For  {110}  oriented  growth,  the 
optimum  initial  growth  temperature  was  970*C.  The 
layer  quality  degraded  rapidly  for  initial  growth  tem¬ 
peratures  lees  than  900*  or  greater  than  580*C.  The 
melt  composition  used  was  1. 4-1.5%  ZnSnPj  with  one 
gram  of  Sn  used  per  run.  Good  growth  was  obtained 
with  cooling  rates  between  2.9*-9*C/hr.  No  layers 
were  grown  on  {110}  oriented  substrates  for  cooling 
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rates  leas  than  U*C/hr.  All  layer*  grown  on  {110} 
oriented  aubetrate*  were  lea*  than  S  *m  thick. 

X-Ray  Diffraction  Measurements 

X-ray  diffraction  measurements  were  made  on  the 
layers  grown  on  {111}  Aa  and  {110}  oriented  GaAa 
using  a  double  crystal  diffractometer  aligned  for  the 
GaAa  333  reflection  for  Cu  Ka  radiation.  Peak  widths 
at  half-peak  intensity  (FWHM)  of  8.7  and  10.6  sec  of 
arc  were  reproducibly  obtained  with  the  {111}  and 
{110}  GaAa  used  for  growth  substrates,  respectively.1 
Some  results  of  these  measurements  are  shown  in 
Fig.  6, 7,  and  8. 

Most  of  the  {111}  As  samples  revealed  a  single 
broadened  GaAs  peak  with  some  slight  background 
reflection  as  seen  in  Fig.  6.  This  suggests  that  this  epi¬ 
taxial  material  lattice  matches  GaAs  to  within  2  X 
10~*A  (4  x  10-%)  based  on  an  8-7  sec  resolution. 
The  FWHM  peak  widths  for  these  layers  were  4S  to 
83  sec.  The  low  intensity  background  reflection  may 
be  due  to  the  rough  interfaces  of  these  layers  or  in¬ 
terfacial  strain. 

Some  of  the  thinner  {111}  As  layers  revealed  a 
secondary  peak  at  about  one  tenth  the  intensity  of 
the  main  peak  separated  by  4.9  min  of  arc  as  seen  in 
Fig.  7.  These  layers  were  3  *m  thick  and  were  grown 
using  a  5’ C/hr  cooling  rate.  Typically,  these  {111}  As 
layers  had  smoother  interfaces  than  those  grown  with 
mailer  cooling  rates,  and  no  broad  background  in¬ 
tensity  was  observed.  The  lattice  constant  associated 
with  the  secondary  peak  in  Fig.  7  is  S.643A  which  is 
0.1%  smaller  than  the  5.65 1A  lattice  constant  ob¬ 
served  by  other  workers  (1,  2)  on  solution- grown 
platelets.  The  FWHM  widths  of  the  main  peaks  were 
20-90  sec,  and  those  of  the  secondary  peaks  were 
around  2  min.  The  observation  of  secondary  peaks 
from  thin  layers  and  broad  background  reflection 
from  thick  layers  suggests  that  there  is  interfacial 
strain  in  these  structures.* 

A  typical  scan  of  the  333  reflection  from  a  layer 
grown  on  {110}  oriented  GaAs  is  shown  in  Fig.  8. 
This  scan  shows  two  approximately  equal  intensity 
peaks,  which  are  separated  by  90  sec.  All  layers  grown 

‘The  SSI  reflection  from  a  (111)  oriented  Ge  analyxor  was  need 
In  aU  nektai  curve*  presented.  This  causes  some  peak  spreading 
due  to  the  different  Lattice  constants  of  the  analyser  sad  sample 
crystals.  The  theoretical  peak  width  for  the  333  redaction  using 
either  Oa  or  GaAs  (111)  oriented  analyser  and  sample  crystals  Is 
4.7  sac  of  are.  IK.  A.  G.  Halllwell,  J.  B.  Childs,  and  S.  O’Hara,  In 
-Gallium  Araenide  and  Related  Compound.,”  C.  BUsum,  Editor,  p. 
SS,  The  Institute  of  Physics,  London  (1973).] 

•The  existence  of  strain  In  these  structures  Is  evidenced  by  the 
broadening  of  the  reflection  peak.  WhUe  strain  In  a  layer  can  ob¬ 
scure  the  rocking  curve  analysis,  the  presence  of  primarily  one 
peak  Indicates  a  lattice  match.  1W.  J.  Bartels  and  H.  Veenvlict,  in 
-Gallium  Araenide  and  Related  Compounds.”  C.  M.  Wolfe,  Editor, 
p.  SIS,  The  Institute  of  Physics,  London  (1979).] 


Rfl.  4.  A  typical  s-rey  dHfructiee  teas  of  the  333  rtf  lectio*  from 
ZaSafj  traps  ea  {111}  A*  arise  ted  GaAs. 


fif.  7.  A  dlffractiaa  scan  of  the  333  rsfiectioa  from  ZnSnF? 
trope  oa  {111}  As  oriented  GaAs  *rowa  with  o  S*C/hr  coolie* 
rate. 


Fig.  I.  A  typical  x-ray  diffraction  scan  from  Zn$nf2  grown  on 
{110}  oriented  GaAs. 

on  {110}  oriented  materiel  revealed  these  two  peaks 
with  the  epitaxial  layer  peak  intensity  scaling  with 
the  layer  thickness.  The  lattice  constant  of  the  layers 
calculated  from  these  data  is  5.6507A  ±  3  x  10~4A 
based  on  a  5.6532A  GaAs  lattice  constant  (4)  which 
is  in  excellent  agreement  with  the  platelet  data. 

Bock  Reflection  Laue  Patterns 

Because  of  the  differences  between  the  diffraction 
measurements  for  {111}  As  and  {110}  oriented  growth, 
back  reflection  Laue  photographs  of  the  layers  were 
made  to  determine  any  differences  in  crystal  struc¬ 
ture.  Typical  photographs  are  shown  in  Fig.  9.  In  all 
photographs,  the  samples  were  aligned  so  that  a 
<110>  axis  of  the  substrate  was  directed  along  the 
incident  x-ray  beam.  A  copper  source  was  used  for 
these  exposures  and  the  sample- to-ftlm  distance  was 
approximately  4  cm. 

Figure  9a  is  a  photograph  of  an  8  ^m  ZnSnPt  layer 
grown  on  {111}  As  oriented  GaAs  whose  diffraction 
scan  was  similar  to  that  shown  in  Fig.  6.  Exposures  of 
this  layer  contain  no  additional  diffraction  spots  com¬ 
pared  with  exposures  of  GaAs.  Hence,  this  ZnSnPj 
layer  has  grown  with  the  sphalerite  structure.  A 
gnomonic  projection  of  the  sphalerite  reflections  ob¬ 
served  in  these  photographs  is  shown  in  Fig.  10a. 
Figure  9b  is  a  photograph  of  a  4  pm  layer  grown  on 
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Fig.  9.  kck  refaction  Loee  photographs  of  ZiSth  layers  grow*  on  (a,  left)  {111}  At  oriented  GaAs  aa 4  (k,  right)  {110}  oriented 
GaAs. 


{110}  oriented  GaAt.  The  diffraction  scan  for  this 
ZnSnPt  layer  is  shown  in  Fig.  8.  The  Laue  photo* 
graph  in  Fig.  9b  shows  a  number  of  diffraction  spots 
which  are  not  present  in  the  GaAs  exposure.  These 
additional  spots  are  due  to  a  higher  degree  of  order¬ 
ing  in  the  ZnSnP2  layer  and  cannot  be  attributed  to 
the  slightly  different  lattice  constant  of  the  layer. 
The  intense  additional  spots  which  are  evident  in 
this  exposure  are  from  high  index  chalcopyrite  167 
and  259  planes,  which  are  listed  in  catalogs  of  powder 
diffraction  patterns.  These  diffraction  spots  are  evi¬ 
dent  despite  their  small  structure  factor  because  they 
meet  the  first-order  Bragg  condition  for  a  wavelength 
of  1.542A,  i.e.,  at  the  high  intensity  Ka  peak  of  the 
Cu  x-ray  spectrum.  These  observations  indicate  that 
the  {110}  oriented  growth  has  the  chalcopyrite  struc¬ 
ture.  A  gnomonic  projection  of  the  chalcopyrite  re¬ 
flections  is  shown  in  Fig.  10b. 


Epitaxial  Orientation 

The  manner  in  which  the  chalcopyrite  structure  of 
the  ZnSnPt  layer  is  oriented  with  respect  to  the 
sphalerite  structure  of  the  GaAs  substrate  is  also  of 
interest  Due  to  the  tetragonal  structure  of  the  epi¬ 
taxial  layer,  there  are  three  ways  in  which  it  may 
orient  relative  to  the  substrate  sphalerite  structure. 
This  is  illustrated  schematically  in  Fig.  11  where 
three  chalcopyrite  conventional  (cubic)  unit  cells  are 
shown  at  three  different  orientations  on  a  {110} 
sphalerite  surface.  In  this  figure,  the  subscript  “c” 
refers  to  chalcopyrite  axes.  The  different  relative 
orientations  are  with  the  [001  ]c  axis  oriented  along 
the  [100],  [010],  or  [001]  axes  of  the  sphalerite  struc¬ 
ture.  Other  researchers  (5)  involved  with  chalcopyrite/ 
sphalerite  epitaxy  have  observed  both  unique  epitaxy 
where  only  one  relative  orientation  is  present  and 
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n«.  II.  Scbowatic  illmtrotioo  of  Hm  rolctWo  oriontation  of  Htf 
chalcopyrito  cmimImmI  unit  colli*  oo  a  {110}  foMorito  sorfoc*. 

mixed  epitaxy  where  more  than  one  relative  orienta¬ 
tion  ia  present. 

Because  ZnSnP*  exhibits  no  tetragonal  distortion, 
the  uniqueness  of  the  epitaxy  cannot  be  determined 
from  the  diffraction  scans  of  the  333/336c  reflections. 
The  chalcopyrite  Laue  spots,  however,  indicate  that 
the  epitaxy  is  mixed  with  the  [001]c  axis  out  of  the 
growth  surface. 

Hall  and  Resistivity  Measurements 

Van  Der  Pauw  resistivity  and  Hall  measurements 
were  made  on  the  {111}  As  and  {110}  layers.  All  layers 
measured  were  found  to  be  p-type.  Layers  grown  on 
{111}  As  oriented  substrates  have  carrier  concentra¬ 
tions  of  1  X  10‘*~1  x  1021  cm-*  and  calculated  mo¬ 
bilities  of  10-70  cmVVsec.  The  resistivity  of  these 
layers  ranged  from  2.5  x  10-4-1.2  x  10-1  Item.  Layers 
grown  on  {110}  oriented  substrates  have  carrier  con¬ 
centrations  of  3  x  10  ‘*-5  x  10‘»  cm-*  with  calculated 
mobilities  of  15-40  cmWsec  and  resistivities  of 
4  x  lO-s-2  x  10-*  ncm.  No  strong  correlations  were 
observed  between  these  data  and  the  growth  con¬ 
ditions  of  the  layers. 

Surface  Effects 

From  the  results  discussed  above,  it  is  evident  that 
the  growth  surface  plays  an  important  role  in  deter¬ 
mining  the  crystal  structure  of  ZnSnPj  grown  on 
GaAs  substrates.  A  similar  observation  was  made  by 
Kroemer  et  al.  (6)  for  GaAs  grown  on  Ge  substrates. 
For  GaAs  on  Ge,  the  {110}  orientation  was  found  to 
be  the  optimum  growth  orientation  for  high  quality 
GaAs  layers.  Kroemer  et  at.  postulated  that  the  {110} 
Ge  surface  is  superior  because  of  the  formation  of  well 
structured  As-like  and  Ga-like  sites.  They  based 
this  conclusion  on  Harrison's  (7)  model  for  the  re¬ 
construction  of  this  nonpolar  surface.  Such  a  recon¬ 
structed  Ge  surface  should  suppress  antiphase  bound¬ 
ary  formation  in  GaAs  producing  layers  with  superior 
properties. 

An  analogous  effect  may  be  responsible  for  the 
growth  of  the  two  ZnSnP*  phases  on  different  GaAs 
substrate  orientations:  The  reconstructed  {110}  GaAs 
surfaces  in  liquid  phase  epitaxy  may  contain  Ga  sites 
which  appear  to  be  Zn-like  or  Sn-like,  while  the  {111} 
As  surfaces  may  not.  In  this  manner,  the  {110}  sur¬ 
faces  would  promote  chalcopyrite  growth  with  anti¬ 
phase  boundaries.  For  instance,  the  {110}  surfaces 
might  reconstruct  into  linear  chalcopyrite-like  sur¬ 
face  chains  in  the  [001]  direction  which  are  out  of 
phase  with  neighboring  chains.  This  would  promote 


the  chalcopyrite  phase  with  antiphase  boundaries,  and 
could  explain  the  somewhat  degraded  electrical  prop¬ 
erties  of  our  chalcopyrite  layers  in  comparison  to  the 
sphalerite  layers.  If  this  is  the  case  then  growing  on 
misoriented  substrates  could  increase  the  lateral 
growth  velocity  and  inhibit  antiphase  boundary  for¬ 
mation.  Another  possibility  is  that  misoriented  sub¬ 
strates  may  produce  a  surface  which  reconstructs  to 
give  longer  range  chalcopyrite-like  regions  which  are 
free  of  antiphase  boundaries. 

Although  Harrison  (7)  considered  the  reconstruc¬ 
tion  of  polar  surfaces  such  as  GaAs,  he  did  not  sepa¬ 
rate  the  metallic  or  nonmetallic  sublattices  into  fur¬ 
ther  sublattices  of  dissimilar  sites  since  such  con¬ 
siderations  were  secondary  to  the  scope  of  his  work. 
Considering  the  experimental  results  discussed  in  this 
paper,  a  study  of  such  reconstruction  may  be  war¬ 
ranted. 

Conclusions 

We  have  developed  a  procedure  for  growing  ZnSnP* 
layers  on  GaAs  in  an  open  tube  liquid  phase  epitaxial 
system.  After  a  preliminary  study  of  ZnSnPg  growth 
on  various  GaAs  substrate  orientations,  the  growth 
procedures  were  refined  for  growth  on  {111}  As  and 
{110}  oriented  substrates.  ZnSnP2  grown  on  {111}  As 
oriented  substrates  with  cooling  rates  smaller  than 
5*  C/hr  was  sphalerite  and  lattice  matched  GaAs. 
ZnSnPg  grown  on  {110}  oriented  substrates  with  cool¬ 
ing  rates  between  2.5* -5*  C/hr  was  mixed  orientation 
chalcopyrite  with  a  lattice  constant  of  5.6507A,  in 
agreement  with  data  on  solution-grown  platelets.  All 
layers  grown  were  p-type  with  carrier  concentrations 
of  1  x  10l*-l  x  10*1  cm-*  and  mobilities  of  10-70 
cmVVsec.  Antiphase  boundaries,  which  are  present 
in  such  mixed  orientation  epitaxy,  may  be  responsible 
for  the  inferior  electrical  characteristics  of  the  chal¬ 
copyrite  layers. 
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Using  back  reflection  Laue  measurements  we  have  examined  the  substrate  orientation  dependence  of  antiphase  domain  boundaries 
for  ZnSnP2  ( c  —  2a)  grown  on  GaAs  by  liquid  phase  epitaxy.  These  results  together  with  those  of  other  workers  indicate  that  (211 } 
substrates  are  best  for  the  growth  of  chalcopyrite  layers  with  c  -  2a,  while  (100)  substrates  are  best  for  those  chalcopyrites  with 
c  <  2a.  The  suppression  of  antiphase  domain  boundaries  can  be  related  to  improved  electrical  properties  for  the  chalcopyrite  layers. 
In  general,  crystal  orientations  with  large  c-axis  components  parallel  to  the  growth  plane  tend  to  be  suppressed.  This  suggests  that 
differences  in  lateral  growth  rate  from  nuclei  with  different  orientations,  rather  than  site  preference,  is  a  dominant  factor  in  the 
suppression  of  antiphase  domain  boundaries. 


1.  Introduction 

The  II-IV-V2  chalcopyrites  have  a  number  of 
interesting  properties  [1)  with  potential  device  ap¬ 
plications.  In  particular,  the  development  of  tech¬ 
niques  by  which  these  materials  can  be  epitaxially 
grown  on  III— V  sphalerite  or  IV  diamond  sub¬ 
strates  can  result  in  heterojunctions  which  com¬ 
bine  nonlinear  optical  properties  with  light  emis¬ 
sion,  detection,  or  signal  processing  capabilities. 
One  of  the  problems  in  achieving  these  heterojunc¬ 
tions,  however,  is  that  the  chalcopyrite  crystal 
structure  is  tetragonal,  whereas  the  sphalerite  and 
diamond  structures  are  cubic. 

The  chalcopyrite  structure  (space  group  I42d)  is 
tetragonal  because  the  II  and  IV  cation  sites  are 
non-equivalent,  while  the  corresponding  sites  in 
the  cubic  sphalerite  and  diamond  structures  are 
equivalent.  Because  of  the  lower  symmetry  of  the 
tetragonal  structure,  many  of  the  chalcopyrites 
exhibit  a  shortened  c-axis  (e/a  <  2  for  the  tetrago¬ 
nal  superlattice)  while  others  do  not.  Conse¬ 
quently.  in  chalcopyrite-on-sphalerite  or  chal- 
copyrite-on-diamond  epitaxy,  there  is  no  unique 

•  Current  address:  Vanan  Associates.  611  Hansen  Way.  Palo 
Alto.  Calfomia  94303.  USA. 


direction  for  the  c-axis,  and  antiphase  domain 
boundaries  are  produced  between  regions  with 
different  c-axis  orientations.  Through  their  strains 
and  contribution  to  scattering  these  antiphase  do¬ 
main  boundaries  are  expected  to  degrade  the 
transport  properties  and  stability  of  the  material, 
and  to  present  problems  in  device  fabrication. 

The  problem  of  antiphase  domain  boundary 
formation  has  received  some  attention  in  the  re¬ 
cent  literature.  Kroemer  et  al.  [2]  have  examined 
the  conditions  under  which  it  should  be  possible  to 
achieve  sphalerite-on-diamond  epitaxy  without  an¬ 
tiphase  domain  boundary  formation.  Considering 
the  results  of  Harrison  et  al.  [3],  they  point  out 
that  there  are  good  reasons  to  expect  substrate 
orientation  to  be  important  in  the  suppression  of 
antiphase  domain  boundaries.  For  example,  the 
{110}  diamond  surfaces  could  coherently  recon¬ 
struct  in  such  a  way  as  to  produce  cation-like  and 
anion-like  sites. 

Wright  et  al.  [4j  found  that  antiphase  domain 
boundaries  are  suppressed  when  sphalerite  GaP  is 
grown  on  a  {211}  surface  of  Si.  They  explain  this 
result  by  demonstrating  that  the  {211}  diamond 
surface  can  consist  of  microscopic  steps  which 
exhibit  cation-like  and  anion-like  sites  without  re¬ 
construction.  These  results  are  also  directly  appli- 
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cable  to  the  reduction  of  antiphase  domain 
boundaries  associated  with  cation-anion  site  inter¬ 
change  in  chalcopyrite-on-diamond  epitaxy.  They 
are  not  applicable,  however,  to  the  suppression  of 
antiphase  domain  boundaries  associated  with  ca¬ 
tion-cation  site  interchange  which  produces  re¬ 
gions  with  different  c-axis  orientations. 

Andrews  et  al.  [5]  have  examined  antiphase 
domain  formation  for  ZnSiAs2  epitaxially  grown 
on  {100}  and  {111}  surfaces  of  GaAs  and  Ge.  In 
this  work  they  were  able  to  determine  growth 
conditions  under  which  the  chalcopyrite  c-axis  grew 
uniquely  normal  to  the  {100}  surfaces.  Thus,  anti¬ 
phase  domain  boundaries  were  suppressed  for  this 
orientation.  For  nominal  {111}  substrates,  how¬ 
ever,  their  results  indicated  that,  of  the  three  possi¬ 
ble  c-axis  orientations,  the  preferred  orientation 
was  the  c-axis  closest  to  the  surface  normal. 

Since  ZnSiAs2  is  one  of  the  II-IV-V2  chal- 
copyrites  which  exhibits  a  foreshortened  c-axis,  the 
question  arises  as  to  whether  it  is  foreshortening  or 
the  tetragonal  atomic  arrangement  which  produces 
the  unique  c-axis  orientation  on  {100}  substrates. 
As  Andrews  et  al.  [5]  point  out,  to  answer  this 
question  it  is  necessary  to  investigate  antiphase 
domain  formation  in  a  chalcopyrite  material  with 
c/a  =  2  such  as  ZnSnP2. 

In  the  present  paper  we  examine  antiphase  do¬ 
mains  in  ZnSnP2  grown  on  {100},  {110},  misori- 
ented  {110},  {111},  and  {211}  surfaces  of  GaAs. 
This  heterojunction  system  is  well  suited  for  such  a 
study  for  the  following  reasons:  (1)  The  GaAs 
substrate  should  suppress  antiphase  domain 
boundaries  associated  with  cation-anion  inter¬ 
change,  leaving  only  those  associated  with  non¬ 
equivalent  cations.  (2)  As  indicated  above,  there  is 
no  foreshortening  in  ZnSnP2  (c  =■*  2a)  so  that  dif¬ 
ferent  substrate  orientations  can  be  examined  with 
relative  ease.  (3)  The  small  lattice  mismatch  of 
0.44%  between  ZnSnP2  (5.6507  A)  and  GaAs 
(5.6532  A)  [6]  minimizes  lattice  strain  and  makes 
epitaxial  growth  easier. 


2.  Growth  procedure 

An  open  tube  liquid  phase  epitaxial  system  was 
used  for  the  growth  of  ZnSnP2  on  GaAs  in  this 


work.  The  details  of  the  system  have  been  previ¬ 
ously  discussed  [6].  The  growth  solutions  were 
composed  of  Zn,  Sn  and  SnP3  synthesized  in  a 
sealed  tube  process  with  a  1.20%  equivalent  ZnSnP2 
concentration  in  Sn.  This  provided  slightly  super¬ 
saturated  solutions  at  the  homogenization  and  ini¬ 
tial  growth  temperature  of  570.0  ±1.0°C,  and 
maintained  melt  saturation  during  a  2  h  homo¬ 
genization  period.  Growth  was  initiated  by  sliding 
the  solution  over  the  substrate  and  cooling  the 
system  linearly  at  3.5  to  4.5°C/h  for  5  to  16  h. 
These  conditions  produced  optimum  quality 
growth  for  all  orientations  examined. 

The  GaAs  substrate  surfaces  were  prepared  as 
follows:  Each  as-received  wafer  was  lapped  and 
chemical-mechanically  polished  using  a  1%  Br- 
CHjOH  solution  until  both  surfaces  were  uniform 
and  free  of  gross  damage.  Such  damage  was  usu¬ 
ally  evident  after  a  20  s  etch  in  the  Br-CH3OH 
solution.  After  scribing  and  breaking  each  wafer 
into  substrates  of  the  proper  dimensions  for  the 
LPE  boat,  the  substrates  were  cleaned  in  boiling 
trichloroethylene,  acetone  and  methyl  alcohol. 
Each  substrate  was  etched  for  2  min  in  a  10  : 1: 1 
H2S04:H202:H20  solution,  10  min  after  initial 
mixing,  and  rinsed  thoroughly  in  deionized  H20 
just  prior  to  loading  the  furnace. 

The  misoriented  {110}  substrates  were  pre¬ 
pared  individually  after  scribing  and  breaking  the 
wafers  into  the  proper  sized  pieces.  Each  {110} 
substrate  was  mounted  on  a  2°  angle  lapping 
block  and  lapped  and  polished  as  above.  The 
substrate  was  then  remounted  face  down  on  a  0° 
lapping  block  and  the  lapping  and  polishing  pro¬ 
cedure  was  repeated.  The  resultant  substrates  had 
parallel  polished  faces  at  an  angle  of  2.0  ±  0.2° 
from  the  initial  {110}  surfaces.  These  substrates 
were  subsequently  cleaned  and  etched  as  indicated 
above. 

The  {110}  substrates  were  misoriented  in  two 
orthogonal  directions:  towards  a  perpendicular 
(110)  and  towards  a  perpendicular  (001)  direc¬ 
tion.  The  direction  of  misorientation  was  de¬ 
termined  from  the  perpendicular  cleavage  planes 
of  the  wafers:  those  misoriented  towards  a  per¬ 
pendicular  cleavage  plane  are  misoriented  towards 
a  perpendicular  (110)  direction;  those  misoriented 
along  a  perpendicular  cleavage  plane  are  misori- 


-9- 


G  A  Davis  el  at.  /  Antiphase  domain  boundary  suppression 


143 


ented  towards  a  perpendicular  (001)  direction. 
For  brevity  these  surfaces  will  be  denoted 
{110}/(110)  and  {110}/(001),  respectively. 


3.  Layer  morphology 

Surface  and  cross-sectional  views  of  a  typical 
ZnSnP2  layer  grown  on  a  {100}  GaAs  substrate 
are  shown  in  fig.  3  of  ref.  [6],  This  layer  and 
subsequent  {100}  layers  had  smooth  interfaces  but 
rough  surfaces  containing  pits  which  often  ex¬ 
tended  to  the  substrate.  Although  we  tried  many 
times,  we  were  not  able  to  improve  the  {100} 
morphology. 

A  typical  {110}  layer  is  shown  in  fig.  5  of  ref. 


Fig.  1.  Surface  and  cross-sectional  views  of  ZnSnP2  grown  on 
GaAs  misoriented  2°  off  the  (110)  towards  a  perpendicular 
<100>. 


[6],  This  ZnSnP2  layer  and  other  {110}  layers 
grown  later  had  smooth  interfaces,  surfaces  with 
some  structure,  and  tended  to  be  somewhat  non- 
uniform  in  thickness. 

Typical  morphology  for  {111}  layers  is  shown 
in  fig.  4  of  ref.  [6].  This  layer  and  subsequent 
{111}  layers  had  very  rough  interfaces  and  smooth 
surfaces.  Back  reflection  Laue  photographs  dem¬ 
onstrated  that  these  layers  were  sphalerite  [6]. 
However,  we  found,  subsequently,  from  electron 
microprobe  analysis  that  all  of  our  {111}  layers 
were  (GaP)2x(ZnSnP2),  _x  alloys,  with  the  Ga  sup¬ 
plied  by  dissolution  of  the  substrate.  Although 
these  are  interesting  results,  they  are  beyond  the 
scope  of  the  present  work. 

Surface  and  cross-sectional  views  of  typical 


Fig.  2.  Surface  and  cross-sectional  view*  of  ZnSnP,  grown  on 
GaAs  muon  ented  2°  off  the  (110)  towards  a  perpendicular 
<110>. 
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Fig  3  Surface  and  cross-sectional  stews  of  ZnSnP.  grown  on 
I  ’1 1 1  As  GaAs 


Fig.  4  Surface  and  cross-sectional  views  of  ZnSnP.  on  |  211  )Ga 
GaAs. 


mtsorienied  {110},.  <001)  ZnSnP.  are  shown  in 
fig.  1.  These  layers  are  5  to  7  gm  thick  for  16  h 
growth  times  with  smooth  interfaces  and  mod- 
eratelv  smtxnh  surfaces  exhibiting  some  structure. 
This  morphology  ts  quite  similar  to  that  achieved 
with  {110}  growth  except  for  somewhat  more 
surface  structure. 

Growth  on  {110}  110)  surfaces  showed 

marked  improvement  over  growth  tin  {110}  and 
I  110}  001-  surfaces.  A  typical  { 1 10}/<1 10) 

ZnSnP  layer  is  shown  in  fig.  2.  These  lasers  are  20 
to  25  gm  thick  for  16  h  growth  times  with  smooth 
interfaces  and  surfaces. 

ZnSnP  lasers  grown  on  211  As  and  21 1  ,Ga 

I 


I 


oriented  GaAs  are  shown  in  figs.  3  and  4.  respec¬ 
tively.  Growth  on  the  {211}  As  orientation  is  20  to 
25  /am  thick  for  a  16  h  growth  time  with  smooth 
interfaces  and  fairly  smooth  surfaces  containing 
some  structure.  The  (211  }Ga  growth  is  also  20  to 
25  /am  thick  for  a  16  h  growth  time  with  smooth 
interfaces  and  surfaces  exhibiting  minimal  struc¬ 
ture.  It  is  interesting  to  note  that  the  cross-sec¬ 
tional  view  of  fig.  4  does  not  show  the  rippling 
observed  in  the  other  cross-sectional  view  suggest¬ 
ing  this  growth  has  improved  crystallinity.  While 
{ 21 1  }Ga  growth  is  the  most  sensitive  to  changes  in 
the  growth  conditions,  it  can  be  reproduciblv  ob¬ 
tained. 
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4.  Structural  measurements 

To  examine  the  structure  of  the  epitaxial  ZnSnP2 
and  determine  the  relative  orientations  of  the  te¬ 
tragonal  axis  of  the  chalcopyrite  layer  with  respect 
to  the  cubic  axes  of  the  sphalerite  substrate,  back 
reflection  Laue  photographs  were  taken.  A  de¬ 
termination  of  the  tetragonal  axis  orientation  rela¬ 
tive  to  the  cubic  axes  is  fundamental  to  this  work 
since  the  existence  of  only  one  relative  orientation 
is  a  necessary  (but  not  sufficient)  condition  for  the 
elimination  of  antiphase  boundaries.  The  tech¬ 
nique  used  was  adapted  from  our  previous  ob¬ 
servations  of  allowed  reflections  from  the  chal¬ 
copyrite  structure  for  {110}  growth  (6). 

The  experimental  arrangement  for  the  photo¬ 
graphs  is  illustrated  schematically  in  fig.  S.  In  all 
Laue  photographs  the  ZnSnP2-GaAs  heterostruc¬ 
tures  were  oriented  so  that  the  Cu  X-rays  were 
incident  along  one  of  the  (110)  directions  of  the 
sphalerite  substrate.  If  the  sample  is  oriented  so 
that  this  corresponds  to  the  [021]  chalcopyrite 
direction,  then  four  additional  diffracted  spots  - 
the  259,  167,  259,  and  167  diffractions  -  will  be 


Fig.  5.  Schematic  of  allowed  chalcopyrite  reflections. 


[021] 


cioo] 

Fig.  6.  Laue  photograph  of  ZnSnP2  grown  on  nominal  (110) 
GaAs  with  primary  X-ray  beam  along  surface  normal. 

evident  in  the  photographs.  It  can  be  seen  from 
fig.  5  that  a  tetragonal  axis  in  the  [001]  direction 
will  produce  these  four  additional  spots  on  one 
side  of  the  photograph.  If  the  chalcopyrite  layer  is 
oriented  with  its  tetragonal  axis  along  the  [010] 
axis  of  fig.  5,  then  the  reflections  would  be  ob¬ 
served  on  the  other  side  of  the  photograph.  Thus, 
with  one  such  Laue  photograph  we  observe  whether 
two  of  the  three  possible  tetragonal  axis  orienta¬ 
tions  ( c-axes)  are  present  in  the  ZnSnP2  layer.  The 
third  possible  relative  orientation  of  the  c-axis  can 
be  observed  by  using  a  different  (110)  direction 
for  the  primary  X-ray  beam.  For  the  case  shown  in 
fig.  5,  one  could  observe  chalcopyrite  growth  with 
the  c-axis  in  the  [100]  direction  by  taking  Laue 
photographs  using  primary  X-rays  in  the  [110]  or 
[101]  directions. 

With  this  procedure  back  reflection  Laue  pat¬ 
terns  were  recorded  for  a  large  number  of  ZnSnP, 
layers  grown  on  {100},  {110},  {110}/(110>, 
{ 1 10}/(001 ),  and  {211}  GaAs  substrates.  (As 
indicated  above,  the  {111}  layers  were  not 


Fig.  7.  Laue  photograph  of  ZnSnP2  grown  on  {110}/<110> 
Ga As  with  primary  X-ray  beam  in  |1 10|  direction  2°  off  the 
surface  normal. 


ZnSnP2.)  Fig.  6  is  a  typical  photograph  for  the 
{110}  growth  and  it  reveals  eight  uniform  chal- 
copyrite  spots.  This  photograph  was  taken  with 
the  primary  X-ray  beam  along  the  surface  normal 
and  indicates  that  {110}  growth  contains  equal 
amounts  of  two  c-axis  orientations. 

Fig.  7  is  a  photograph  of  a  ZnSnP,  layer  grown 
on  {110}/(110)  GaAs.  The  direction  of  the 
primary  beam  in  this  case  was  the  [110]  axis  2° 
from  the  surface  normal.  This  photograph  con¬ 
tains  four  high  intensity  chalcopyrite  spots  associ¬ 
ated  with  one  c-axis  orientation  as  well  as  four  low 
intensity  chalcopyrite  spots  (opposite  the  former 
spots)  associated  with  a  less  abundant  second  c-axis 
orientation. 

A  typical  Laue  pattern  for  layers  grown  on 
{211}  substrates  is  shown  in  fig.  8.  For  a  surface 
normal  of  (1 12]  the  primary  X-ray  beam  was  in  the 
(101]  direction.  This  photograph  shows  the  four 
asymmetric  spots  observed  when  only  one  c-axis 


Fig.  8.  Laue  photograph  of  ZnSnP2  grown  on  (112)  GaAs  with 
primary  X-ray  beam  in  f  101 )  direction. 


orientation  is  present  in  the  epitaxial  layer.  The 
absence  of  other  orientations  was  confirmed  by 
photographs  with  the  primary  beam  in  the  other 
(110)  directions.  (In  one  direction  only  sphalerite 
spots  are  observed.) 


[100) 

Fig.  9  Schematic  illustration  of  (he  relative  orientations  of 
ZnSnP,  grown  on  (100)  GaAs. 


G.  A.  Davis  el  a!.  /  Antiphase  domain  boundary  suppression 


147 


[100] 

Fig.  10.  Schematic  illustration  of  the  observed  relative  orienta¬ 
tions  of  ZnSnP2  grown  on  {110}  GaAs. 


The  results  of  these  Laue  measurements  for  all 
of  the  ZnSnP2  layers  are  summarized  schematically 
in  figs.  9  to  12.  In  these  figures  the  shaded  planes 
indicate  the  growth  surfaces  and  the  arrows  along 
the  axes  of  the  cubic  unit  cell  indicate  the  relative 
abundance  of  the  three  possible  chalcopyrite  c-axis 
orientations.  The  arrows  are  denoted  [001  ]c  to 
indicate  the  chalcopyrite  tetragonal  axis.  The 
lengths  of  the  arrows  illustrate  the  relative  abun¬ 
dance  of  each  orientation. 

As  indicated  in  fig.  9,  all  three  possible  c-axis 
orientations  are  observed  in  {100}  growth  with  the 
c-axis  normal  to  the  surface  most  abundant. 

Fig.  10  shows  the  relative  c-axis  orientations 
observed  in  {110}  growth.  In  this  case  the  epi- 


[100] 

Fig.  11.  Schematic  illustration  of  the  relative  orientations  of 
ZnSnPj  grown  on  (110}  GaAs  misoriented  2°  towards  a  per¬ 
pendicular  (110). 


[100] 


Fig.  12.  Schematic  illustration  of  the  relative  orientation  of 
ZnSnP2  grown  on  (211 }  GaAs. 

taxial  layers  contained  equal  amounts  of  material 
with  tetragonal  axes  in  the  [100]  and  [010]  direc¬ 
tions  but  none  in  the  [001]  direction. 

There  was  no  detectable  difference  between  the 
Laue  photographs  for  (110}/(001)  growth  and 
that  shown  in  fig.  6  for  {110}  growth.  Photographs 
of  {110}/<110)  growth,  however,  revealed  a  sig¬ 
nificant  increase  in  the  abundance  of  material  with 
the  c-axis  in  the  [010]  direction  and  a  decrease  in 
the  amount  with  the  c-axis  in  the  [100]  direction. 
This  is  shown  in  fig.  11.  In  some  layers  grown  on 
the  {110}/(110>  surfaces,  the  smaller  relative  c- 
axis  orientation  was  not  observed.  Variations  in 
the  relative  abundance  of  the  two  c-axis  orienta¬ 
tions  may  be  due  to  unintentional  misorientation 
of  the  substrates  during  preparation. 

The  {211}As  and  {211}Ga  growth  revealed  the 
best  results  in  that  only  one  c-axis  orientation  was 
observed.  This  is  illustrated  in  fig.  12  which  shows 
the  observed  relative  orientation  of  the  tetragonal 
axis  along  the  [001]  direction  for  a  (112)  surface. 
This  habit  was  observed  in  all  {211}  growth,  re¬ 
gardless  of  morphology. 


5.  Compositional  analysis 

To  determine  the  composition  of  the  epitaxial 
layers,  electron  microprobe  analysis  was  per¬ 
formed  for  each  substrate  orientation.  With  the 
exception  of  the  { 1 1 1 }  layers,  all  were  stoichiomet¬ 
ric  ZnSnPj  within  the  accuracy  of  the  analysis 


-14- 


14# 


G.A. Davis  el  at.  /  Antiphase  domain  boundary  suppression 


which  was  about  1%.  The  analysis  indicated  a 
slight  tendency  to  incorporate  Ga  from  the  GaAs 
substrate  into  the  layers  near  the  interface  with  the 
Ga  level  being  less  than  1  mol%.  No  As  was  found 
in  any  of  these  layers. 


6.  Electrical 


Most  of  the  ZnSnP2  layers  were  grown  on 
semi-insulating  GaAs  substrates  to  facilitate  elec¬ 
trical  characterization.  Van  der  Pauw  measure¬ 
ments  were  then  performed  on  the  various  orienta¬ 
tions.  except  for  the  (100)  samples  which  had 
poor  morphology,  and  the  carrier  concentrations 
and  mobilities  were  calculated  with  a  simple  one- 
band  model. 

The  layers  grown  on  (110}  and  (110}/(001) 
substrates,  which  were  composed  of  material  with 
two  c-axis  orientations  of  roughly  equal  abun¬ 
dance,  had  hole  concentrations  of  3  to  5  X  1014 


cm'3  and  mobilities  of  15  to  40  cn^/V-s.  The 


layers  grown  on  (110}/(110)  substrates,  with  one 
of  the  two  c-axis  orientations  suppressed,  and  those 
grown  on  {211}  substrates,  with  only  one  c-axis, 
had  hole  concentrations  of  2  to  5  X  1018  cm  * 3  and 
mobilities  of  20  to  45  cn^/V  •  s.  Thus,  the  suppres¬ 
sion  of  antiphase  domain  boundaries  apparently 
yields  improved  electrical  transport  properties. 


7.  Conclusions 


A  comparison  of  our  results  for  ZnSnP2  (c  -  2a ) 
grown  on  GaAs  to  those  of  Andrews  et  al.  (5)  for 
ZnSiAs2  (c<2a)  grown  on  GaAs  indicates  cer¬ 
tain  orientation  preferences  for  the  suppression  of 
antiphase  domain  boundaries  in  chalcopyrite-on- 
sphalerite  epitaxy.  For  the  chalcopynte  com¬ 
pounds  with  c  -  2a,  {211}  substrates  are  ap¬ 
parently  the  best,  while  for  those  with  c  <  2a. 
( 100}  substrates  are  better  for  the  suppression  of 
antiphase  domain  boundaries.  Here,  of  course,  we 
are  referring  to  antiphase  domain  boundaries  asso¬ 
ciated  with  material  having  different  c-axes.  There 
may  still  be  antiphase  domains  which  conserve  the 
c-axis,  such  as  those  which  are  produced  by  rota¬ 
tion  of  90°  about  the  c-axis. 

Our  results  and  those  of  Andrews  et  al.  [5] 


indicate  that  antiphase  domains  with  the  largest 
c-axis  components  in  the  growth  plane  tend  to  be 
suppressed.  The  misoriented  {110}  experiments 
demonstrate  this  rather  convincingly.  From  fig.  10 
it  can  be  seen  that  tilting  a  (110)  surface  towards 
the  perpendicular  [001]  direction  does  not  change 
the  relative  components  of  [100]  and  [010]  c-axis 
orientations  in  the  growth  plane,  while  from  fig.  1 1 
tilting  a  (110)  towards  the  perpendicular  [IlO] 
does. 

Although  we  do  not  understand  the  reasons  for 
the  orientation  dependence  of  antiphase  domain 
boundary  suppression,  it  appears  to  be  more  than 
site  preference.  Slight  surface  misorientations  do 
not  affect  the  results,  so  surface  reconstruction  is 
probably  not  a  factor.  Also,  considering  the  results 
for  {211}  substrates,  there  is  no  obvious  reason 
why  this  orientation  would  distinguish  between 
cation  sites  in  a  manner  which  would  minic  the 
chaicopyrite  structure.  Although  there  may  be  some 
energetically  favorable  site  preference,  it  seems 
more  likely  that,  somehow,  the  high  lateral  growth 
velocity  suppresses  the  formation  of  antiphase  do¬ 
mains. 
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G.A.  DAVIS 

Semiconductor  Research  Laboratory,  Washington  University,  St.  Louis,  MO  63130 


ABSTRACT 

Misordering  defects  which  occur  in  the  growth  of  chalcopyrite/sphalerite 
systems  have  been  investigated  through  the  LPE  growth  of  ZnSnP2  on  GaAs. 

X-ray  diffraction  measurements  reveal  the  presence  of  misordering  defects 
associated  with  multiple  tetragonal  orientations  in  the  epitaxy.  Growth  on 
{ 211 }  oriented  substrates  reveals  the  lowest  concentrations  of  orientation- 
related  defects  while  growth  on  vicinal  {110  }  oriented  substrates  indicates 
the  strongest  cation  site  differentiation  required  to  produce  long-range 
ordering.  The  use  of  critically  misoriented  {110}  substrates  is  proposed 
for  the  growth  of  highly  ordered  chalcopyrite  layers  on  sphalerite  substrates. 


INTRODUCTION 

The  chalcopyrite-phase  semiconducting  compounds  have  generated  much 
interest  in  the  past  for  application  of  their  unique  optical  and  electronic 
transport  properties  [1].  Much  of  this  interest  has  moderated  recently  in 
light  of  the  substantial  difficulties  encountered  in  the  growth  of  large 
high-quality  crystals  suitable  for  device  processing:  solution- grown  and 
vapor-grown  platelets,  from  which  most  of  the  data  concerning  the  chalco- 
pyrites  has  been  determined,  contain  significant  impurity  concentrations 
in  addition  to  being  inconsistent  with  current  planar  processing  techniques. 
Epitaxial  growth  of  the  chalcopyrites  on  their  sphalerite  analogs  has  re¬ 
sulted  in  layers  containing  large  concentrations  of  defects  which  severely 
degrade  the  properties  of  these  materials.  Pending  the  development  of 
procedures  capable  of  producing  low  defect  concentration  epitaxial  layers, 
these  materials  will  remain  technologically  insignificant.  The  objective 
of  this  work  is  to  study  the  nature  of  chalcopyrite/sphalerite  epitaxy  in 
order  to  develop  growth  conditions  which  suppress  defect  formation. 


DEFECT  CONSIDERATIONS 

Defects  are  formed  in  the  growth  of  mixed  epitaxial  systems  when  the 
substrate  exhibits  a  lower  degree  of  ordering  than  that  exhibited  by  the 
epitaxial  material.  As  a  consequence,  the  substrate  can  not  ensure  proper 
long-range  ordering  of  the  epitaxial  constituents  and  misordering  defects 
result. 

This  type  of  defect  was  initially  discussed  by  Holt  [2]  for  the 
sphalerite/diamond-structure  system.  In  this  system  the  defects  are  in 
the  form  of  boundaries  across  which  there  exists  a  cation/anion  site 
interchange  (incoherence)  and  are  called  antiphase  domain  boundaries  (ADBs). 
The  nature  and  suppression  of  ADBs  have  been  extensively  studied  in  the 
GaP/Si,  GaAs/Ge  and  GaAs/Si  systems  [3-6]  for  applications  to  sphalerite/ 
diamond- structure  device  integration.  ADB  suppression  has  been  pursued 
through  the  use  of  diamond-structure  substrates  which  promote  cation/anion 
site  differentiation  consistent  with  the  sphalerite  structure.  Such  site 
differentiation  may  be  encouraged  through  reconstruction  of  the  substrate 
surface  [3,7]  into  dissimilar  cation-like  and  anion-like  sites  or  through 
dissimilar  bond  coordination  of  the  sites  [4].  This  latter  method  has 
resulted  in  MBE  GaP  free  of  ADBs  for  growth  on  {211  }  oriented  Si  substrates. 
Single  domain  OMVPE  GaAs  grown  on  {  100  }  oriented  Si  substrates  demonstrates 


the  additional  facility  of  highly  constraining  growth  conditions  in  elim¬ 
inating  ADBs.  Sphalerite-phase  ordering  is  not  supported  by  reconstruction 
or  bond  coordination  for  the  {100}  orientation  and  there  are  additional 
complications  associated  with  charge  accumulation  and  dipole  formation  at 
such  a  polar  interface  (8J.  It  is  noted  that  these  electrostatic  complica¬ 
tions  may  be  relieved  through  exchange  of  anions  and  cations  across  the 
heterointerface . 

The  situation  one  encounters  in  chalcopyrite/sphalerite  growth  is 
somewhat  more  complicated  due  to  the  number  of  constituents  involved  and 
the  complexity  of  the  cation  sublattices.  Presumably  there  will  be  no 
defects  generated  from  cation/anion  ambiguity  since  the  well  ordered  anion 
sublattice  of  the  substrate  is  consistent  with  that  of  the  epitaxial  mater¬ 
ial.  Hence  we  will  only  be  concerned  with  defects  arising  from  misordering 
of  the  cation  sublattices.  We  may  distinguish  between  two  types  of  defects 
in  this  system;  (1)  antiphase  domain  boundaries  being  those  defects  caused 
by  a  simple  translation  of  the  epitaxial  structure  which  results  in  cation 
incoherence,  and  (2)  orientation  domain  boundaries  (ODBs)  being  those 
defects  at  the  boundaries  between  regions  having  different  tetragonal 
orientations.  These  are  illustrated  in  Figure  1  for  hypothetical  growth 
on  {100}  oriented  sphalerite  substrates.  ODBs  are  easily  inferred  from 


(a)  (b) 


Figure  1  (a)  Antiphase  domain  boundary  and  (b)  Orientation  domain  boundary 
for  hypothetical  {100}  orientation  growth.  The  conventional  cell 
has  been  shifted  in  the  right  side  of  (a)  for  clarity. 

simple  diffraction  measurements  and  their  existence  indicates  ambiguous 
or  inconsistent  cation  ordering  at  the  heterointerface.  Since  different 
tetragonal  orientations  present  different  atomic  structures  at  the  hetero¬ 
interface,  the  predominance  of  one  tetragonal  orientation  indicates 
preferential  cation  coordination  which  is  at  least  short-ranged.  The 
approach  to  eliminating  these  defects  is  identical  to  that  used  to  eliminate 
ADBs  in  the  sphalerite/diamond-structure  systems:  by  ensuring  proper 
cation  site  differentiation  at  the  heterointerface  both  ADBs  and  ODBs  will 
be  eliminated.  The  problems  encountered  are  also  identical  and  one  must  be 
cautious  to  select  non-polar  (with  respect  to  the  cations)  surfaces  for 
growth  or  ensure  adequate  redistribution  at  the  heterointerface  to  alleviate 
charge  accumulation  and  dipole  formation.  It  is  noted  that  the  bond  coordi¬ 
nation  approach  can  not  be  utilized  in  this  system  (for  planar  growth 
surfaces)  since  all  substrate  cations  are  identically  coordinated. 


EXPERIMENTAL 


The  LPE  growth  and  X-ray  diffraction  characterization  of  ZnsnP2  on 
GaAs  [9,10]  was  used  for  these  defect  suppression  studies.  ZnSnP2  is  of 
only  minor  interest  for  device  applications  due  to  the  small  cation 
electronegativity  difference  and  consequent  weak  anisotropy  it  exhibits. 

It  is  of  primary  interest  to  this  study,  however,  due  to  this  lack  of 
cation  distinction  and  proclivity  to  disorder. 

The  LPE  growth  procedure  was  implemented  in  an  open  tube  system  util¬ 
izing  a  standard  horizontal  sliding  boat.  The  growth  solutions  were 
comprised  of  Sn,  Zn,  and  SnP3  with  an  equivalent  composition  of  1.20Z 
ZnSnP2  in  Sn  providing  a  saturated  solution  at  the  570°C  initial  growth 
temperature.  Cooling  rates  of  3. 5-4. 5  C/hr  for  5-16  hrs  produced  epitaxial 
layers  4-30um  thick  depending  upon  substrate  orientation.  The  growth 
process  was  self-limiting  for  those  orientations  yielding  thin  layers, 
notably  the  {100},  {110},  and  {111}  orientations,  due  to  spontaneous 
nucleation  in  the  melt.  {100},  {110},  misoriented  {110},  {111},  and  {211} 
substrate  orientations  were  used  in  this  study.  Photomicrographs  of  a 
superior  epitaxial  layer  are  shown  in  Figure  2. 


(a)  (b) 


Figure  2  (a)  Surface  and  (b)  Oilineated  cross-sectional  views  of  a 

ZnSnP2  layer  grown  on  a  misoriented  {110}  GaAs  substrate. 

Back-reflection  Laue  photographs  were  used  to  determine  the  relative 
abundance  of  the  three  possible  tetragonal  orientations  and  hence  the 
existence  of  ODBs.  In  these  photographs  the  samples  were  oriented  such 
that  the  0.5 ram  Cu-target  primary  X-ray  beam  was  directed  along  a  <110> 
direction  with  a  sample-to-f ilm  distance  of  4  cm  (35Z  of  the  long  dimension 
of  the  photographs).  Epitaxial  material  having  a  <021>  direction  along 
the  primary  beam  produces  four  asymmetric  reflections  in  addition  to  the 
sphalerite  reflections.  These  reflections,  the  259,  259,  167,  and  167, 
are  easily  observed  despite  their  small  structure  factors  since  they  meet 
the  first  order  Bragg  condition  for  the  Ka  peak  of  the  Cu  X-ray  spectrum. 

A  typical  photograph  revealing  one  relative  orientation  of  the  tetragonal 
axis  is  shown  in  Figure  3.  Through  photographing  two  orthogonal  <110> 
directions,  the  relative  abundance  of  the  three  possible  tetragonal 
orientations  may  be  determined. 
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Figure  3  Laue  photograph  of  ZnSnP2  grown  on  a  misoriented  {110}GaAs 
substrate.  The  four  asymmetric  high- intensity  refractions, 
allowed  for  the  chalcopyrite  structure,  indicate  a  single 
tetragonal  orientation  having  the  [021]  direction  normal 
to  the  surface. 


DISCUSSION 

The  observed  abundance  of  the  tetragonal  orientations  is  illustrated 
in  Figure  4  for  growth  on  {100},  {211},  {110},  and  misoriented  {110}  sub¬ 
strates.  In  this  figure  the  shaded  planes  indicate  the  substrate  surfaces 
and  the  arrows  along  the  conventional  axes  of  the  substrate  lattice  indi¬ 
cate  the  relative  abundance  of  material  with  the  tetragonal  axis  oriented 
along  the  indicated  direction.  Epitaxial  layers  grown  on  {111}  oriented 
substrates  are  sphalerite-phase  alloys  of  ZnSnP2  and  GaP  attributed  to 
an  instability  of  this  orientation  under  the  growth  conditions.  The  {100} 
orientation  growth  exhibits  all  possible  tetragonal  orientations  with  the 
surface  normal  orientation  being  slightly  more  abundant,  as  illustrated  in 
Figure  4  (a).  This  growth  habit  suggests  that  cation  site  differentiation 
is  sufficiently  weak  so  that  no  tetragonal  orientation  dominates  and  hence 
no  long-range  ordering  is  supported.  In  terms  of  suppressing  ODBs,  the  use 
of  {211}  oriented  substrates  yields  the  best  results  with  epitaxial  layers 
consistently  exhibiting  only  one  tetragonal  orientation,  as  illustrated  in 
Figure  4  (b).  While  this  growth  habit  indicates  preferential  localized 
cation  coordination,  it  is  not  clear  whether  ordered  cation  distinction 
extends  across  the  substrate  surface  as  is  necessary  to  eliminate  ADBs  and 
produce  single  domain  growth.  Layers  grown  on  {110}  oriented  substrates 
reveal  only  two  of  the  tetragonal  orientations,  with  a  complete  absence  of 
material  oriented  with  its  tetragonal  axis  in  the  surface  plane,  as  illus¬ 
trated  in  Figure  4  (c).  The  equal  abundance  of  the  two  out-of-surface 
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Figure  4  Observed  abundance  of  tetragonal  orientations  of  (a)  {100}, 
(b)  {112},  (c)  {110},  and  (d)  misoriented  {110}  substrates. 
The  shaded  planes  indicate  the  substrate  surfaces. 


orientations  is  observed  in  all  Laue  photographs  and  suggests  that  the 
single  crystalline  grains  are  significantly  smaller  than  the  area  of  the 
primary  X-ray  beam.  The  prohibitively  large  strain  and  electrostatic 
energies  associated  with  such  a  large  concentration  of  ODBs  (>>2xl0^cm2) 
suggest  that  the  cation  site  differentiation  of  this  orientation  is  sub¬ 
stantial.  Subsequent  experiments  utilizing  substrates  misoriented  2°  off 
of  a  {110}  towards  a  perpendicular  <110>  reveal  a  reduction  of  one  of  the 
two  tetragonal  orientations  observed  in  nominal  {110}  growth,  as  illustrated 
in  Figure  4  (d).  In  these  layers  the  in-surface  orientation  is  not  observed. 
Also,  some  Laue  photographs  reveal  only  one  tetragonal  orientation,  indi¬ 
cating  an  increase  in  the  grain  size.  Through  the  use  of  substrates  criti¬ 
cally  misoriented  in  this  manner,  it  may  be  possible  to  induce  long-range 
chalcopyrite-phase  ordering. 

One  possible  explanation  of  the  {110}  orientation  results  is  the  ex¬ 
istence  of  a  marked  preference  for  one  cation  species  to  occupy  all 
available  cation  sites  at  the  heterointerface.  Such  a  site  preference  is 
somewhat  compatible  with  the  two  out-of-surface  orientations  which  contain 
alternating  large  II-V  and  IV-V  regions,  and  incompatible  with  the  in¬ 
surface  orientation  which  contains  alternating  cations  at  the  surface.  One 
would  anticipate  such  a  preference  of  one  species  over  another  based  on  the 
different  bond  strengths  of  these  cations  at  the  heterointerface.  This 
behavior  is  apparently  observed  in  the  ordering  of  AlxGai-xAs  into  a  mono- 
layer  superlattice  structure  when  grown  on  {110}  oriented  substrates  [11], 
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wherein  the  cation  electronegativity  difference  is  extremely  small.  Such 
a  planar  cation  preference  would  produce  defect-free  growth  on  {021}  ori¬ 
ented  substrates  as  illustrated  in  Figure  5.  The  corresponding  {011}  chal- 
copyrite  orientation  exhibits  alternating  II-V  and  IV-V  layers  not  exhibited 
by  the  other  corresponding  orientations.  As  with  {100}  sphalerite/diamond- 
structure  growth,  this  configuration  is  electrostatically  prohibitive  unless 
there  is  an  exchange  of  constituents  across  the  hetero inter face. 
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Figure  5  (Oil)  Oriented  chalcopyrite  structure  on  (021)  oriented 
sphalerite  substrate  illustrating  the  layered  nature  of 
this  orientation.  The  conventional  cells  of  the  substrate 
and  epitaxy  are  denoted  by  fine  dashed  lines. 

Finally,  it  is  interesting  to  note  the  correlation  between  the  degree 
of  disordering  and  the  self- limiting  thichnesses  of  the  epitaxial  layers. 
Those  orientations  which  reveal  a  high  degree  of  disordering,  the  {100}, 
{110},  and  {111},  are  typically  on  the  order  of  5pm  thick  while  those  which 
reveal  a  low  degree  of  disordering,  the  {211}  and  misoriented  {110}  are 
significantly  thicker,  being  on  the  order  of  30pm  thick.  This  behavior 
indicates  a  misordering-related  impediment  to  the  continued  growth  at  the 
sample  surface  due,  for  instance,  to  electrostatic  contributions  which 
prohibit  the  incorporation  of  cations  in  the  regions  of  misorientation 
defects. 


CONCLUSIONS 

The  suppression  of  misordering  defects  has  been  studied  through  the 
LPE  growth  of  ZnSnP2  on  GaAs.  While  growth  on  the  {211}  orientations 
exhibits  the  lowest  concentration  of  ODBs,  the  use  of  vicinal  {110}  sub¬ 
strates  appears  to  be  the  most  promising  for  the  elimination  of  both  ADBs 
and  ODBs  due  to  their  strong  cation  differentiation.  In  particular,  the 
use  of  {021}  oriented  substrates  may  promote  long-range  ordering  based  on 
the  results  of  vicinal  {110}  growth  experiments.  A  further  investigation 
on  the  concentration  and  form  of  misordering  defects  through  the  use  of 
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imaging  diffraction  techniques  would  clearly  be  useful  to  complete 
this  study.  The  application  of  constraining  growth  conditions,  avail¬ 
able  in  growth  procedures  other  than  LPE,  may  be  of  additional  utility. 
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Abstract — The  mixing  enthalpy  of  ternary  tetrahedral  semiconductor  alloys  is  fairly  well  described  by 
regular  solution  theory,  with  a  thermodynamic  interaction  parameter  that  is  sensitive  to  the  lattice  spacing 
of  the  binary  constituents.  We  derive  an  estimate  of  the  interaction  parameter  from  a  model  which  ascribes 
the  mixing  enthalpy  to  bond  distortions  associated  with  the  alloy  formation,  and  relates  these  to  the 
macroscopic  elastic  properties  of  the  crystal.  Numerical  estimates  are  given  for  the  18  alloys  with  cations 
Al,  Ga,  In  and  anions  P,  As,  Sb  and  are  compared  with  experimental  values  and  alternative  models.  To 
within  a  single  adjustable  parameter,  the  predictions  agree  with  experiment  and  are  consistent  with  those 
of  the  delta  lattice  parameter  (DLP)  model.  A  further  calculation  of  the  elastic  energy  associated  with 
composition  fluctuations  (clustering)  in  these  alloys  indicates  that  this  energy  is  suffident  to  suppress 
clustering  above  the  critical  mixing  temperature. 


1.  INTRODUCTION 

The  calculation  of  phase  diagrams  and  of  other 
thermodynamic  properties  of  tetrahedral  ternary 
III-V  semiconductor  alloys  of  the  form  AXB,_XC  (or 
AC,DX  .,)  using  regular  solution  theory  has  proved  to 
be  a  useful  and  reasonably  accurate  approximation. 
Several  of  the  ternary  alloys  with  cations  Al,  Ga,  In 
and  anions  P,  As,  Sb  have  important  present  or 
potential  applications  in  electronic  technology.  These 
eighteen  alloys  are  mixtures  of  binary  compounds 
with,  in  general,  different  lattice  parameters.  The 
lattice  parameters  of  the  binary  compounds  are 
different  when  the  AC  and  BC  bond  lengths  are 
different,  and  the  lattice  parameter  of  the  alloy  then 
differs  from  both.  The  lattice  parameters  of  all  the 
ternary  alloys  that  can  be  formed  from  the  ions  listed 
except  those  involving  the  cation  pair  Al,  Ga  are 
sensitive  to  composition.  One  then  may  take  the  view 
that  the  bonds  are  distorted  or  rearranged  in  forming 
the  alloy.  The  energy  associated  with  the  bond  dis¬ 
tortion  can  be  estimated  from  the  macroscopic  elastic 
properties  of  the  crystals.  We  discuss  the  contribution 
of  this  energy  to  the  enthalpy  of  mixing  and  its  effect 
on  clustering  and  spinodal  decomposition.  We  sug¬ 
gest  that  it  is  a  major  or  dominant  factor  in  both 
phenomena. 

2.  MIXING  ENTHALPY 

In  regular  solution  theory  the  molar  enthalpy  of 
mixing  of  a  binary  solution  is  given  by 


AHm  =  x(  1  -x)G 


where  G  is  the  interaction  parameter  and  where  x  and 
I  -  x  are  the  fractions  of  the  components  making  up 
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the  solution.  In  terms  of  a  picture  of  pairwise  inter¬ 
action  energies,  the  interaction  parameter  is 


G  =  6N& 


(2) 


where  6 N0  is  the  number  of  interacting  pairs  on  the 
fee  sublattice  and 
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(3) 


(D 


is  the  excess  of  the  AA  and  BB  bond  energies  between 
like  ions  over  the  AB  bond  energy  of  unlike  ions.  In 
the  ternary  alloys  under  discussion  these  second 
nearest  neighbor  interactions  are  thought  to  be  rather 
small,  consistent  with  the  small  experimental  values 
of  the  interaction  parameter. 

Stringfellow,  using  the  covalent  bonding  ideas  of 
Phillips  and  Van  Vechten[l],  has  introduced  the  delta 
lattice  parameter  model  [2]  which  provides  values  of 
the  interaction  parameter  of  all  the  alloys  in  terms  of 
a  single  adjustable  parameter.  The  agreement  be¬ 
tween  the  DLP  theory  and  experiment  is  remarkably 
good  for  a  theory  with  a  single  adjustable  parameter. 
The  model  takes  a  less  specific  approach  than  the 
identification  of  eqns  (2)  and  (3)  of  the  interaction 
parameter  with  the  sum  of  binary  bond  energies,  and 
assumes  that  the  energy  is  a  more  global  property  of 
the  semiconductor  that  depends  only  on  the  lattice 
parameter  of  the  crystal.  The  idea  is  that  the  bonding 
energy  is  proportional  to  the  band  gap  which,  in  turn, 
is  proportional  to  a  power  of  the  lattice  parameter. 
This  is  similar  to  but  distinct  from  Van  Vechtens 
proposal  [3]  that  G  is  produced  by  the  reduction  of  the 
band  gaps  due  to  the  disorder  associated  with 
alloying.  The  DLP  model  leads  to  an  approximate 
value  for  Si  of 


G,  =  (35/8)  KA2J~ 


(4) 
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where  A  «  a  —  b,  J  —  1/2 (a  4-  6),  with  a  and  b  the 
lattice  parameters  of  -4C  and  BC  respectively.  The 
quantity  K  is  obtained  by  a  least  squares  fit  to  the 
experimental  data  and  is  equal  to  1.15  x  I07cal/mole 
A25.  The  excellent  fit  to  the  experimental  data  de¬ 
pends  largely  on  the  A 2  dependence  of  Q. 

We  wish  to  show  here  that  if  part  or  all  of  Q  is 
interpreted  as  being  associated  with  strain  or  lattice 
parameter  mismatch,  then  it  is  readily  shown  that  this 
contribution  exhibits  the  same  d2  dependence. 

It  is  observed  that  the  average  lattice  spacing  of  the 
alloys  closely  follows  Vegard's  law 

4)  =  xa  +(1  -x)b.  (5) 

Since  the  difference  between  the  lattice  parameters  a 
and  b  of  the  binary  compounds  corresponds  to  the 
different  AC  and  BC  bond  lengths,  the  bonds  must  be 
distorted  or  rearranged  to  form  the  alloy.  If  a  and  b 
do  not  differ  too  widely,  the  energy  of  the  rearranged 
configuration  must  be  of  second  order  in  a  -  d,  and 
b  —1^.  Thus  the  contribution  to  the  enthalpy  of 
formation  of  the  alloy  is 


AHm 


X 

7 


(xf^-a^  +  O  -xXdt-bf), 


(6) 


where  X  is  an  adjustable  parameter.  With  d,  from  eqn 


(5)  this  becomes 

AHm  -  jc<  I  -x)X  (7) 

with  the  same  A1  dependence  as  the  DLP  model. 

Since  we  are  viewing  the  effect  as  an  elastic  dis¬ 
tortion,  we  can  obtain  a  crude  numerical  estimate  for 
X  from  the  macroscopic  elastic  properties  of  the 
crystal.  At  one  extreme  one  can  assume  a  strict  virtual 
crystal  model  of  the  alloy,  which  would  identify  the 
bond  distortion  energy  with  that  of  on  isotropic 
compression  or  dilatation.  At  x  »  1/2,  the  principal 
strains  are  given  by  l/2d  and  the  elastic  energy  per 
unit  volume  is 

(8) 

yielding 

o-.i (|J  ,9, 

where  B  is  the  bulk  modulus  and  V  is  the  molar 
volume  of  the  alloy  with  x  —  1/2.  Table  I  lists  Q,  and 
O'  for  the  ternary  compounds  under  discussion. 

Equation  (9),  which  contains  no  adjustable  param¬ 
eters,  yields  higher  interaction  parameter  values  than 


Table  I.  Experimental  and  calculated  values  of  the  interaction  parameter.  Bulk  moduli  by  interpolation 
between  binary  compound  values  from  Refs.  [8]  and  |9fc  interaction  parameters  (in  k  cal/mole)  in  columns 
3. 4,  and  S  quoted  from  Stringfcllow,  Ref.  (2)  (small  differences  in  the  calculated  values  in  column  4  because 
we  use  fewer  significant  numbers  than  Ref.  [2]  for  A  it) 
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the  DLP  theory.  A  much  better  fit  to  the  experimental 
data  is  obtained  if  the  prediction  of  eqn  (9)  is  reduced 
by  a  factor  of  about  4.4,  resulting  in  the  column 
labeled  Q'  in  Table  I.  Indeed,  the  entries  in  that 
column  fit  the  observed  values  as  closely  as  String* 
fellow's[2]  “best  fit”  2.45  power  law.  Moreover  a 
plausible  case  can  be  argued  for  a  reduction  of  the 
estimate  by  a  factor  of  this  order. 

It  has  been  demonstrated  recently  [4]  that  although 
the  average  lattice  parameter  closely  follows  Vegard’s 
law,  the  AC  and  BC  bond  lengths  do  not  take  the 
common  nearest  neighbor  value  of  the  virtual  crystal 
but  remain  much  closer,  by  a  factor  of  4  or  5,  to  their 
lengths  in  the  binary  compounds.  If  it  is  assumed  that 
the  stretching  and  compression  of  nearest  neighbor 
bonds  is  the  only  distortion,  this  reduces  the  estimate 
of  >1  by  the  square  of  this  factor,  of  order  20.  This, 
however,  almost  certainly  underestimates  the  effect  of 
the  distortion,  since  the  conservation  of  nearest 
neighbor  bond  lengths  requires  a  distortion  of  the 
bond  angles  and  of  second  neighbor  bond  lengths 
and  angles. 

In  addition,  the  fact  that  the  AC  and  BC  bond 
lengths  tend  to  remain  constant  (to  about  80%) 
independent  of  alloy  composition,  necessarily  implies 
that  the  A  and  B  atoms  in  the  alloy  are  much  more 
correlated  than  in  a  random  solution.  These  strong 
correlations  could  well  be  due  to  the  strong  elastic 
forces.  Neither  our  theory  nor  any  other  theory 
currently  takes  this  into  account  properly.  Thus 
within  the  scope  of  our  model.  Table  1  furnishes 
strong  evidence  that  this  quasi-elastic  energy  is  a 
significant,  and  perhaps  a  dominant  part  of  the 
enthalpy  of  mixing. 

X  CLUSTERING 

A  binary  mixture  whose  enthalpy  of  mixing  is 
positive  has  a  tendency  to  separate  into  its  com¬ 
ponents.  This  tendency  manifests  itself  in  the  form  of 
spinodal  decomposition  and  the  development  of  a 
miscibility  gap  below  a  critical  mixing  temperature 
and  in  the  form  of  clustering  on  a  microscopic  scale 
above  this  temperature.  In  the  alloys  under  discussion 
here,  this  tendency  is  most  readily  interpreted  in 
terms  of  preferential  pairwise  interactions  as  ex¬ 
pressed  by  eqns  (2)  and  (3),  but  it  does  not  depend 
on  this  interpretation.  Jones,  Porod,  and  Ferry  [5] 
have  recently  examined  the  degree  of  clustering  that 
may  be  expected  to  occur  above  the  critical  tem¬ 
perature  of  mixing  as  a  function  of  the  interaction 
parameter  interpreted  in  this  manner,  of  the  tem¬ 
perature,  and  of  the  mole  fraction  jc.  We  extend  this 
calculation  to  include  an  elastic  energy  which  is  large 
enough  to  overcome  the  clustering  tendency. 

Since  the  lattice  parameters  of  the  alloys  follow 
Vegard’s  law,  the  ions  of  the  alloy  at  least  on  the 
average  occupy  positions  on  a  virtual  crystal  lattice 
with  the  lattice  parameter  of  eqn  (5).  Note  that  this 
statement  is  not  in  conflict  with  the  observation  cited 
above(4]  that  the  bond  lengths  need  not  follow 
Vegard’s  law. 
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Gustering  produces  a  fluctuation  from  the  average 
composition,  while  the  region  of  the  cluster  must 
occupy  the  volume  prescribed  by  the  average  virtual 
crystal  lattice.  Thus  the  region  of  the  cluster  is 
strained.  The  energy  associated  with  the  strain  tends 
to  reduce  the  degree  of  clustering,  but  we  shall  see 
that  it  does  not  eliminate  a  miscibility  gap  below  the 
critical  mixing  temperature.  The  effect  is  similar  to 
the  “lattice  latching”  [6]  and  stabilization  against 
decomposition [7]  associated  with  the  mismatch 
strains  of  heteroepitaxy.  In  the  present  instance  we 
deal  with  the  strains  arising  from  fluctuations  in  a 
homogeneous  medium.  In  either  case,  the  existence  of 
a  strained  region  presupposes  the  existence  of  the 
crystal  lattice  to  be  matched  by  the  continuing 
growth.  In  epitaxy,  this  lattice  is  furnished  by  the 
substrate.  In  homogeneous  growth  it  depends  on  the 
formation  of  a  nucleus,  which  at  a  given  composition 
can  only  form  at  temperatures  above  the  spinodal 
decomposition  curve.  This  is  the  reason  why  the 
strain  energy  affects  clustering  but  not  the  miscibility 
gap. 

Let  the  average  composition  of  an  alloy  be  A^B^ 

C  and  the  local  composition  AZB,_XC.  Then  the 
lattice  parameter  is 

<4  =  *.«+(  l-x0)b.  (10) 

If  the  crystal  were  permitted  to  relax  everywhere,  the 
local  lattice  parameter  would  be 

d  =  xa  +(1  -  x)b.  (11) 

However,  the  compatibility  constraint  imposes  a  lo¬ 
cal  strain 

6  =rf-4,  =  (*-x0Xa-&)s{d,  (12) 

the  second  equality  defining  the  composition 
fluctuation  {. 

The  energy  density  of  the  strain  can  be  computed 
as  in  eqn  (8)  and  is 


per  unit  volume. 

Unlike  eqn  (8)  for  the  elastic  contribution  to  the 
mixing  enthalpy,  eqn  (12)  should  not  be  a  significant 
overestimate  of  the  strain  energy.  In  the  present 
application  we  are  computing  the  energy  of  actual 
compressions  and  dilatations  in  terms  of  measured 
elastic  constants  and  the  nature  of  the  bond  dis¬ 
tortions  accompanying  the  strain  is  irrelevant.  There¬ 
fore,  if  the  mixing  enthalpy  is  dominated  by  the  bond 
distortion  effects  discussed  in  Section  2,  then  it  is 
evident  that  the  strain  energy  of  composition 
fluctuations  is  more  than  adequate  to  overcome  the 
clustering  tendency. 

A  similar  argument  can  be  made  if  the  clustering 
is  ascribed  to  preferential  pairwise  interactions.  The 
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molar  bond  energy  density  of  the  ternary  alloy, 
taking  into  account  the  4  nearest  and  12  next  nearest 
neighbor  interactions  is 

£  *  A/'0(4£^*  +4£j(1  —  x) 

I2t^»x(l -x)  +  6€m(1  -x)1).  (14) 

where  iA  and  t,  are  the  AC  and  BC  interaction 
energies.  The  contribution  to  the  bond  energy  density 
of  a  composition  fluctuation  {  is 
£(*<>  +  £)  -  E(x o)  -  n.  given  by 

7  -  m*A  ~  +  >2<4.  -  1 2* „  -  1  2x„£  K  -  *) 

(15) 

with  £  as  defined  in  eqn  (3). 

In  a  crystal  of  a  given  composition,  the  con¬ 
tribution  to  the  interaction  energy  linear  in  the 
fluctuations  must  by  definition  average  to  zero. 
Therefore  only  the  term  in  {J  need  be  considered 
from  the  purpose  of  minimizing  the  free  energy.  But 
from  eqn  (2)  this  term  is  just 

tf  =  -Qt\  (16) 

so  that  the  net  molar  energy  of  the  fluctuations,  from 
eqns  (13)  and  (16),  is 

WV  +  r,  =(Q'-QKJ  (17) 

averaged  over  the  crystal. 

Equation  (17)  yields  a  positive  (or  zero)  fluctuation 
energy  for  each  of  the  ternary  alloys  we  are  consid¬ 
ering,  indicating  that  clustering  is  suppressed  by  the 
strain  energy. 

4.  CONCLUSION 

The  view  is  represented  in  our  model  that  the 
formation  of  a  tetrahedral  semiconductor  alloy  is 


accompanied  by  bond  distortion,  and  that  an  esti¬ 
mate  of  the  electronic  energy  of  this  distortion  can  be 
obtained  from  the  macroscopic  elastic  properties  of 
the  crystal.  Superficially  this  view  appears  distinct 
from  the  basis  of  the  DLP  model,  which  is  deduced 
from  the  empirical  lattice  parameter  dependence  of 
the  optical  band  gap  and  from  Van  Vechtens  pro¬ 
posal.  But  the  bandgap  itself,  as  pointed  out  by 
Phillips  and  Van  Vechten[l],  is  a  measure  of  the 
covalent  bonding  energy.  The  elastic  constants  of  a 
crystal  depend  as  directly  on  the  bonding  energy  as 
its  optical  and  dielectric  properties.  Therefore  it  is 
perhaps  not  too  surprising  that  models  derived  from 
optical  and  from  elastic  crystal  parameters  should 
yield  predictions  of  similar  form  for  the  thermo¬ 
dynamic  properties.  The  merit  we  may  claim  for  the 
present  model  is  that  it  contributes  a  new  and  con¬ 
venient  experimental  approach  to  the  thermodynamic 
quantities  and  that  it  may  yield  a  fruitful  theoretical 
insight. 
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Spatial  composition  fluctuations  in  semiconductor  alloys  4xfl,_xC  involve  elastic  enerty  when  the  lattice 
parameters  of  AC  and  BC  differ.  Therefore,  the  fluctuations  are  expected  to  be  reduced  below  those  ex¬ 
hibited  by  a  purely  random  distribution  of  the  constituent  ions.  We  estimate  the  extent  of  the  reduction, 
and  relate  it  to  recent  observations  of  Raman  linewidths  and  line  shapes  in  Gaxln|_xAs  and  Ga^A^.^As. 


Ternary  tetrahedrally  coordinated  semiconductor  alloys 
AXB\-XC  or  ACXD\-X,  where  A  and  B  are  group  III  (or  II) 
elements,  and  C  and  D  group  V  (or  VI)  elements  can  for 
many  purposes  be  described  as  binary  mixtures  whose  com¬ 
ponent  particles  reside  (approximately)  on  the  sites  of  an  fee 
lattice.  Specifically,  thermodynamic  parameters  such  as 
mixing  enthalpy1  and  spinodal  curves2  conform  fairly  well  to 
the  rules  describing  regular  solutions. 

The  key  quantity  of  regular  solution  theory  is  the  interac¬ 
tion  parameter  fl ,  which  can  be  interpreted  in  terms  of  pair 
interaction  energies;  thus 

ft  ”  -jNeZdEAt  —  Ea*  —  Ebb)  . 

where  No  is  Avogadro’s  number,  Z  the  coordination 
number  (12  for  fee),  and  EtJ  are  the  (negative)  pair  interac¬ 
tion  energies.  A  negative  value  of  fl  indicates  a  tendency 
to  ordering  (compound  formation),  a  positive  value  a  ten¬ 
dency  to  decomposition.  All  the  observed  interaction 
parameters  of  the  alloy  group  with  A,B  —  Al.Ga.In; 
(Zn,Cd,Hg)  C.D  -  P,As,Sb;  (Se.Te)  are  positive  or  zero. 
Thus,  one  expects — and  observes — miscibility  gaps  and 
spinodal  decomposition  below  a  critical  mixing  temperature 
characteristic  of  each  alloy. 

The  observed  interaction  parameters  are  found  to  be  pro¬ 
portional  to  the  square  of  the  lattice  parameter  difference 
between  the  pure  compounds  that  constitute  the  alloy.  This 
dependence  has  been  interpreted,  in  the  so-called  delta- 
lattice-parameter  (DLP)  model  in  terms  of  the  lattice 
parameter  dependence  of  the  average  band  gap,2  and  more 
recently,  especially  after  details  of  the  composition  depen¬ 
dence  of  bond  lengths  in  Gai_zInxAs  were  revealed  by  ex¬ 
tended  x-ray  absorption  fine  structure  (EXAFS)  measure¬ 
ments,4  in  terms  of  the  bond  distortions  associated  with  the 
formation  of  the  alloy.  By  relating  the  bond  distortion  ener¬ 
gies  to  the  macroscopic  elastic  constants  of  the  crystals,  a 
good  case  could  be  made  for  ascribing  most  or  all  of  the  ob¬ 
served  mixing  enthalpy  to  this  cause.5,6 

A  positive  mixing  enthalpy  would  normally  produce  clus¬ 
tering  of  like  particles  above  the  critical  mixing  temperature, 
and  an  estimate  of  the  degree  of  clustering  to  be  expected 
was  published  recently.7  For  the  model  that  was  used  in 
this  calculation,  the  origin  of  the  positive  interaction  param¬ 
eter  is  irrelevant.  It  has  been  pointed  out,  however,  that  in 
the  type  of  “lattice  gas"  that  represents  these  alloys,  with  a 
composition-dependent  lattice  parameter,  if  an  alloy  crystal 
is  formed  at  all,  elastic  interactions  will  suppress  the  cluster¬ 
ing.5  It  was  also  suggested  that  these  elastic  interactions 
were  likely  to  produce  an  ordering  tendency,  that  is  to  say, 


that  local  composition  fluctuations  might  be  reduced  below 
the  level  to  be  expected  of  a  purely  random  distribution  of 
the  constituents.  What  follows  is  a  simple  thermodynamic 
model  for  the  composition  correlations  induced  by  the  elas¬ 
tic  energy. 

The  lattice  parameters  of  these  alloys  follow  Vegard’s  law 
quite  accurately.  Any  deviation  from  a  linear  dependence  of 
lattice  parameter  on  composition  (bowing)  is  small  enough 
to  be  quite  negligible  for  the  purpose  of  this  argument. 
Consider  a  region  of  the  alloy  large  enough  to  be  regarded 
as  macroscopic.  If  the  composition  of  this  region  deviates 
from  the  average  composition  of  the  alloy,  but  must  be  ac¬ 
commodated  in  the  volume  of  the  virtual  lattice  that  “be¬ 
longs  to  it,”  it  will  suffer  compression  or  dilatation. 

The  energy  associated  with  the  strain  is  readily  comput¬ 
ed.5  Let  the  average  composition  of  the  alloy  be  AXJ3\-XJC 
and  the  local  composition  AXB\-XC.  Then  the  lattice 

parameter  is 

do**xga  +(l-x0)6  ,  (1) 

where  a  is  the  lattice  parameter  of  AC  and  b  the  lattice 
parameter  of  BC.  If  the  crystal  were  permitted  to  relax 
everywhere,  the  local  lattice  parameter  would  be 

d~xa+(\-x)b  \  (2) 

however,  the  compatibility  constraint  imposes  a  local  strain 

8-rf-do—  (x-xo)(a  -b)mAx(a  —  b)  ,  (3) 

with  an  associated  energy  density  (per  unit  volume) 

»'(Ax)j-18b|-2-^|  (Ax)J  ,  (4) 

[a  +b  J 

where  B  is  the  bulk  modulus  of  the  crystal. 

Equation  (4)  is  based  on  isotropic  compression  or  dila¬ 
tion;  an  analogous  expression  can  be  obtained  for  a  planar 
strain.  In  that  form  the  elastic  energy  is  thought  to  be 
responsible  for  the  “lattice  tatching”  observed  in  heteroepi- 
taxy,*  the  growth  of  a  single  lattice-matched  alloy  composi¬ 
tion  over  a  range  of  melt  compositions.  It  also  accounts  for 
the  successful  heteroepitaxial  growth  of  alloy  compositions 
unstable  in  the  bulk  at  the  growth  temperature. 

It  has  been  shown5  that  a  reasonable  fit  to  the  interaction 
parameters  of  the  III-V  alloys  considered  here  is  given  by 
the  semiempirical  formula 

IK'- 0.226  x  18b|-j-^|J  (5) 
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making  a  contribution  to  the  mixing  enthalpy  per  mole  of 

A//-  W‘x( \-x)V  .  (6) 

where  V  is  the  molar  volume  of  the  alloy. 

The  entropy  of  mixing  is  approximately 


S  -  In 


Na\N,\ 


NAA\N„\{NAt\Y 


(7) 


where  the  N,  are  numbers  of  ions  and  the  Nu  numbers  of 
second  neighbor  pairs.  Per  mole  this  is 

S  =•  /M  x  lax  +  ( 1  -  x  )ln(  1  -  x )  -  x2  ln(x2) 

—  (1— x)2ln(l-x)2— 2x(l— x)lnx(l— x)l  .  (8) 


Using  Eqs.  (4),  (6),  and  (8)  in  the  thermodynamic  formula’ 
for  macroscopic  fluctuations  of  concentration  on  an  assem¬ 
bly  of  N  sites  with  total  Gibbs  free  energy  G 


(  (Ajx  )2) 

we  find 


kT 

(dJ<//a*J)  r.r.N 


(9) 


<(A,)2)-ix(1-x)[l4-^^>-^(l-x) 


(10a) 

or,  in  terms  of  the  critical  mixing  temperature  Te  given  by 
IRTc—  W'V  —  0.226  WV 

<(Aor)J)-  — j|+137(J./T)jf(1_jc)  j  .  (10b) 

One  may  guess  conservatively  that  the  disorder  grown 
into  a  crystal  corresponds  roughly  to  a  value  frozen  in  at  the 
growth  temperature.  Thus,  if  we  compare  a  lattice-matched 
alloy  such  as  Gao.jAl«.jAs  with  Tt  =*  0  with  the  mismatched 
alloy  In&jGaojAs  with  Te  —  600  K  grown  epitaxially  at 
T  —  900  K,  we  expect,  from  Eq.  (10a)  or  (10b)  to  find  a  re¬ 
lative  reduction  of  the  mean-square  concentration  fluctua¬ 
tions  in  the  mismatched  alloy  by  a  factor  —  0.3.  The  reduc¬ 
tion  would  be  expected  to  be  more  significant,  depending  on 
cooling  rate  and  on  the  growth  parameters,  if  appreciable 
solid-state  diffusion  takes  place  after  growth. 

Some  qualitative  experimental  support  for  this  prediction 
is  to  be  found  in  linewidths  and  line  shapes  observed  by 
Parayanthal  and  Poliak10  in  Raman  scattering  from 
AI,rGai-xAs  and  InxGai-xAs.  In  this  work  it  was  found 
that  the  broadening  and  asymmetry  of  the  Raman  line  was 
more  pronounced  in  the  lattice-matched  alloy  AlxGai_xAs 
than  in  the  mismatched  InxGai-xAs. 

A  “spatial  correlation"  model  has  been  successfully  used 
to  interpret  the  effect  of  microcrystalline,11  and  implantation 
damage12  disorder  on  Raman  linewidths  and  line  shapes. 
This  model  replaces  the  phonons  of  the  infinite  ordered 
crystal  with  phonons  localized  within  “correlation  regions” 
whose  size  is  a  measure  of  the  partial  ordering  in  the  imper¬ 
fect  crystal.  As  a  result,  the  q  -  0  momentum  selection  rule 
of  the  perfect  crystal  is  relaxed,  and  phonons  in  a  region  of 
the  BriUouin  zone  corresponding  to  the  correlation  length 
participate  in  the  scattering.  Experimental  values  of  the 
correlation  length  we  deduced  by  fitting  the  observed 


linewidths  and  line  shapes,  assuming  a  Gaussian  size  distri¬ 
bution  of  the  correlation  regions,  generally  in  reasonable 
agreement  with  the  estimated  sizes  of  crystallites11  or  un¬ 
damaged  regions.12 

Substitutional  disorder  in  an  alloy  similarly  breaks  the 
crystal’s  translational  symmetry  and,  hence,  the  q  —  0  selec¬ 
tion  rule  for  Raman  scattering.  Since  the  perturbation  is  the 
phonon  modulation  of  the  dielectric  tensor,  the  effective 
size  of  the  scattering  regions  is  given  by  the  correlation 
length  of  this  tensor's  spatial  fluctuations.11  No  theoretical 
prediction  has  so  far  been  made  for  this  quantity,  but  Paray¬ 
anthal  and  Pollack10  were  able  to  fit  their  observed  Raman 
lines  with  a  spatial  correlation  model  in  which  the  correla¬ 
tion  length  is  simply  a  phenomenological  parameter  deter¬ 
mined  from  the  experimental  data.  For  their  epitaxially 
grown  alloys,  this  parameter  ranges  from  about  8  to  above 
30  lattice  parameters. 

Of  the  two  alloy  series  for  which  results  are  given, 
GaxAI|-IAs  has,  within  experimental  accuracy,  zero  mixing 
enthalpy,  and  the  cations  should  constitute  a  perfect  lattice 
gas  on  their  fee  sublattice.  Their  concentration  fluctuations, 
which  may  be  taken  as  measure  of  the  disorder,  then  are 
given  by  Eqs.  (10)  with  7V-0.  The  Raman  lines  bear  this 
out  qualitatively,  in  the  sense  that  the  linewidths  and  asym¬ 
metries  are  greatest  for  samples  with  x  new  O.S. 

The  broadening  of  the  Raman  lines  even  in  the  complete¬ 
ly  random  alloy  GaxAl|-xAs  is  modest  as  compared  with 
the  spectrum  of  amorphous  materials,  and  the  spatial  corre¬ 
lation  model  is  successful  in  the  interpretation  of  the  line 
shapes.  This  suggests  that  one  might  assign  the  role  of  the 
defects  that  disrupt  the  phonon  modes  to  fairly  large  fluc¬ 
tuations  from  the  average  composition.  The  correlation  or 
coherence  length  can  then  be  related  to  the  statistics  of  large 
fluctuations. 

In  a  random  binary  alloy,  the  probability  distribution  of 
ions  is  binomial.  With  a  coordination  number  Z  -  12,  the 
distribution  of  even  the  nearest-neighbor  shell  is  fairly  well 
represented  by  a  Gausaian,  and  the  probability  of  oc¬ 
currence  of  a  deviation  5s  (  is  approximated  by 

/*({)  — erfc(i)  .  (11) 

where  (■f/[2((Ajt)2)  11/2. 

In  order  to  arrive  at  an  order  of  magnitude  estimate  of 
the  correlation  length,  we  take  the  size  of  the  defect  to 
comprise  one  cation  shell  of  Z  ions.  If  a  deviation  of  rela¬ 
tive  magnitude  f  is  to  occur  in  a  region  containing  NZ  ca¬ 
tions,  we  should  have  NZP{()~  1>  and,  therefore,  the 
average  distance  between  such  defects,  which  we  take  to  be 
the  correlation  length,  is 

L-Nula-[ZtTfc(())~ula  .  (12) 

where  a  is  the  lattice  parameter.  Equation  (12)  yields 
values  consistent  with  the  estimates  based  on  the  Raman 
spectra  if  {  is  about  two  standwd  deviations.  Using  the 
average  value  of  L  for  the  GaAs-like  mode  in  GaxAI|.xAs 
for  the  four  samples  with  x  new  O.S,  we  find  f  -  1.88. 

Equations  (10)  and  (12)  also  provide  a  comparison  of  the 
correlation  lengths  in  GaxIni_xAs  and  GaxAI|-xAs.  Be¬ 
cause  of  the  suppression  of  concentration  fluctuations,  we 
expect  the  scattering  centers  associated  with  wide  fluctua¬ 
tions  to  be  rarer,  increasing  the  correlation  length.  Using 
the  factor  0.3  estimated  from  Eq.  (10)  for  the  reduction  in 
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the  mean  square  fluctuation  in  Eq.  (12),  we  find 

- 1; .  in) 

£°S.sAVsA* 

This  factor  is  larger  than  the  numbers  deduced  from  the  Ra¬ 
man  spectra  by  Parayanthal  and  Poliak,  suggesting  that  addi¬ 
tional  phonon  scattering  processes  are  active. 

The  reduction  in  randomness  discussed  here  is  expected 


to  have  other  experimentally  testable  consequences,  which 
will  be  considered  elsewhere. 

I  am  indebted  to  Dr.  A.  Sher  for  suggesting  this  study  and 
for  several  illuminating  discussions,  to  Dr.  F.  Poliak  for 
prepublication  information  about  his  results,  and  to  Dr.  C. 
M.  Wolfe  for  many  helpful  comments.  This  research  was 
sponsored  by  the  Air  Force  Office  of  Scientific  Research, 
Air  Force  Systems  Command,  (J.  S.  Air  Force  under  Grant 
No.  AFOSR82-0231. 
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7 .  Mow  band-gap  photoresponse  of  In,  _  ,Gay  P-GaAs  heterojunctions 
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Semiconductor  Research  Laboratory,  Washington  University,  Sl  Louis,  Missouri  63130 
(Received  9  July  1984;  accepted  for  publication  7  September  1984) 

Below  band-gap  photovoltaic  response  with  a  half-power  point  at  1.31  eV  is  observed  in  p- type 
In,  ^Oa,P  on  n-type  GaAs  heterojunctions  with  { 1 1 1 )  interfaces.  This  response  is  apparently 
due  to  photon-assisted  tunneling  of  carriers  across  the  interfacial  energy  gap,  which  is  smaller 
than  the  energy  gaps  of  the  constituent  materials.  The  same  mechanism  can  be  employed  in  other 
heterojunctions  with  appropriate  energy-band  lineups  to  obtain  longer  wavelength  infrared 
response. 


Several  schemes  have  been  used  in  intrinsic  infrared 
photodetectors  to  obtain  response  below  the  band  gap  of  the 
constituent  semiconductors.  One  of  these  methods  is  to  form 
an  alloy  between  two  compounds  which  have  conduction 


and  valence  bands  of  opposite  symmetry,  such  as  PbTe  and 
SnTe.'  Since  the  conduction  and  valence  bands  must  inter¬ 
change  at  some  alloy  composition,  they  move  toward  one 
another  producing  alloy  band  gaps  smaller  than  the  end- 
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points.  Another  method  uses  electroabsorption  (photon-as¬ 
sisted  tunneling)  to  shift  the  absorption  edge  to  longer  wave 
lengths.2  This  has  been  employed  in  reverse-biased  GaAs' 
and  Ga,  ,  Al,  As  devices.4  A  recently  proposed  scheme  is  to 
form  an  alloy  between  two  materials  with  different  crystal 
structures.  In  the  region  where  the  alloy  undergoes  a  phase 
change  from  one  structure  to  another,  the  absorption  edge  is 
lowered  presumably  due  to  antisite  disorder.  This  behavior 
has  been  observed  in  (GaAs),  ,Ge2,  alloys.'  In  the  present 
work  we  discuss  and  present  experimental  evidence  for  what 
appears  to  be  another  method  of  obtaining  longer  wave¬ 
length  response. 

The  scheme  is  essentially  the  inverse  of  a  photon  emis¬ 
sion  process  previously  described.'1  It  employs  the  electric 
Held  at  the  interface  of  a  heterojunction  to  tunnel  charge 
carriers  across  an  interfacial  energy  gap  smaller  than  the 
energy  gaps  of  the  semiconductors  on  either  side.  Photore¬ 
sponse  is  obtained  by  the  photoexcitation  of  electrons  from 
the  valence  band  of  one  material  to  virtual  states,  and  then 
electron  tunneling  into  the  conduction  band  of  the  other  ma¬ 
terial.  Conversely,  electrons  from  the  valence  band  of  one 
material  can  tunnel  into  virtual  states  and  then  be  photoex- 
cited  into  the  conduction  band  of  the  other  material.  The 
latter  process  is  equivalent  to  hole  tunneling,  and  both  are 
illustrated  schematically  in  Fig.  1 .  The  formalism  for  these 
physical  processes  is  that  of  band-to-band  electroabsorption 
in  a  nonuniform  electric  field  which  has  been  previously  ex¬ 
amined  in  some  detail.  '  4'7 

One  of  the  heterojunction  systems  which  exhibits  the 
required  energy  band  lineup  is  In,  .Ga, P-GaAs.  An 
abrupt  heteroj unction  model  for  this  system  is  shown  in  Fig. 
1.  Several  of  the  parameters  of  the  model  vary  from  sample 
to  sample  and  were  determined  in  the  following  ways:  the 
energy  gap  of  GaAs  was  measured  by  photoconductivity  on 
high-purity  epitaxial  samples  to  be  1.39  cV.  The  energy  gap 
of  In,  ,Ga,  P  was  determined  from  its  dependence  on  lat¬ 
tice  constant,1  where  the  latter  was  obtained  by  double-crys¬ 
tal  x-ray  diffraction  scans.  On  this  basis  the  energy  gap  of  the 
lattice-matched  alloy  Ino  ^Gso  ,,  P  was  1. 90  eV.  The  separa¬ 
tions  between  the  electrochemical  potential  energy,  and 
the  band  edges  were  determined  from  Hall  measurements: 


FIG.  I.  Abrupt  heterojunction  model  for  In,  ,Ga,  P-GaAs  indicating 
schematically  the  photon-assisted  tunneling  processes  for  electrons  on  the  p 
side  and  holes  on  the  n  side 
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first  on  the  GaAs  substrates  and  then  on  In,  .Ga,  P  layers 
grown  on  semi-insulating  GaAs  by  liquid  phase  epitaxy  un¬ 
der  the  same  conditions  as  those  grown  on  n  ’  substrates. 
The  built-in  potential  and  the  band  bending  on  either  side 
were  obtained  from  capacitance-voltage  measurements. 
From  these  measurements  all  the  parameters  of  the  model, 
including  Ac, ,  Je,  .  and  At,,  were  determined  experimen¬ 
tally 

For  a  lattice-matched  heterojunction  with  a  nominal 
1 1 1 1 1  interface,  the  band  offsets  were  Atv  =  0.08  eV,  Acc 
=  0.39  eV,  and  Ac,  =  1.31  eV.  Although  these  results  are 
consistently  obtained  experimentally,  it  should  be  pointed 
out  that  they  are  not  close  to  the  expected  theoretical  values. 
Using  a  dielectric  model,4  for  example,  we  calculate  ioniza¬ 
tion  energies  for  GaAs  and  lno4«Gao ,,  P  of  5.68  eV  and  5.91 
eV,  respectively.  These  ionization  energies  and  the  experi¬ 
mental  energy  gaps  yield  band  offsets  of  Acv  =  —  0.23  eV, 
Aec  =  0.28  eV,  and  Ac,  equals  the  band  gap  of  GaAs.  Thus, 
below  band-gap  photoresponse  is  not  predicted  from  this 
dielectric  model.  Although  there  is  probably  no  reason  to 
expect  good  agreement  between  theory  and  experiment,  we 
have  some  preliminary  evidence  that  the  valence-band  offset 
is  in  the  predicted  direction  for  In,  .,Ga,  P-GaAs  hetero¬ 
junctions  with  nominally  charge-neutral 10  1 2 1 1 1  interfaces. 
Thus,  at  least  part  of  the  difference  in  band  lineups  may  be 
due  to  dipoles  associated  with  the  residual  charge  at  (111) 
interfaces. 

Figure  2  shows  the  room-temperature  photovoltaic  re¬ 
sponse  for  a  1 1 1 1 1  heterojunctibu  at  0  and  —  0.5  V  applied 
bias.  (These  diodes  were  rather  leaky  and  exhibited  substan¬ 
tially  increased  noise  at  higher  reverse  voltages.)  For  these 
measurements  the  device  was  biased  with  a  constant-voltage 
power  supply  and  light  was  incident  on  a  cleaved  surface 
perpendicular  to  the  heterojunction  interface.  The  salient 
features  are  a  below  band-gap  peak  with  a  half-power  point 
at  1.31  eV  which  increases  substantially  with  reverse  bias  as 
expected  for  tunneling.  The  1.31-eV  half-power  point  is 
equal  to  the  value  obtained  for  the  interfacial  energy  gap, 
Ac,.  This  exact  agreement,  however,  should  not  be  taken 
seriously.  Also,  in  the  heterojunction  the  In,  _  .Ga.Pepitax- 
ial  layer  was  doped  to  produce  2.5  x  1017  cm  ~ '  holes,  while 
the  GaAs  substrate  has  1.4  x  1011  cm  “  ’  electrons.  Thus, 


4W0T0M  ■HOT  («V) 


FIG.  2  Responstvity  a*  a  function  of  photon  energy  for  the  heterojunction 
of  Fig.  I  showing  below  band-gap  response  at  two  bias  voltages 
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most  of  the  potential  was  in  the  In,  ,  Ga,  P  where  it  should 
enhance  the  lower  effective  mass  and  higher  probability  elec¬ 
tron  tunneling.  The  GaAs  absorption  edge  is  apparently 
shifted  to  higher  than  band-gap  energy  due  to  the  high  elec¬ 


tron  concentration. 


To  test  the  concept  that  the  below  band-gap  response  is 
associated  with  tunneling  across  the  interfacial  energy  gap, 
several  other  experiments  were  performed.  When  the  doping 
of  the  p- type  In,  _  ,Ga,P  was  increased  from  2.5  x  10'7  to 
6.0x  10'*  cm-3  on  the  same  n-type  GaAs  substrates,  the 
jelow  band-gap  response  was  greatly  reduced.  Increased 
doping  on  thep  side  is  expected  to  reduce  electron  tunneling. 
Since  large  effective  mass  hole  tunneling  from  the  n  side  is 
much  less  probable,  the  smaller  below  band-gap  response  is 
consistent  with  a  tunneling  mechanism.  To  eliminate  the 
possibility  that  the  below  band-gap  response  was  due  to  im¬ 
purity  transitions  in  the  space-charge  region  on  either  side  of 
the  interface,  photoconductivity  measurements  were  per¬ 
formed  on  the  epitaxial  layer  and  on  the  substrate.  These 
measurements  indicated  no  below  band-gap  photoresponse 
at  the  same  energy  due  to  impurities.  Except  for  the  coinci¬ 
dence  of  the  response  energy  with  Ae,,  we  cannot,  however, 
rule  out  the  possibility  of  impurity  transitions  at  the  inter¬ 
face  itself. 

In  conclusion,  we  have  observed  below  band-gap  pho¬ 
tovoltaic  response  in  In,  _  ,  Ga,  P-GaAs p-n  heterojunctions 
which  is  apparently  due  to  photon-assisted  tunneling  across 
the  interfacial  energy  gap  between  the  two  materials.  Al¬ 
though  the  observed  response  for  |  111 )  interfaces  at  1.31  eV 


is  not  at  a  particularly  useful  wavelength,  it  should  be  possi¬ 
ble  to  apply  the  same  idea  for  infrared  detectors  at  longer 
wavelengths.  We  are  currently  investigating  the  possibility 
of  obtaining  smaller  interfacial  energy  gaps  in  the 
In,  _ ,  Ga,  P-GaAs  system  by  examining  interfaces  with  dif¬ 
ferent  orientations.  Other  lattice-matched  heterojunction 
systems,  of  course,  should  also  be  examined.  Our  results  in¬ 
dicate,  however,  that  heterojunctions  with  the  required  band 
lineups  cannot  be  predicted  with  great  reliability. 

This  research  was  sponsored  by  the  Air  Force  Office  of 
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a.  Below  Bandgap  Emission  and  Absorption  in 
ZnSnP2/GaAs  Heterojunctions 
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I.  INTRODUCTION 

TPHERE  IS  CONSIDERABLE  current  interest  in  HI-V 
compound  and  alloy  beceroscractnres  because  of  tbe 
capability  they  provide  for  tailoring  energy  bands  in 
lattice-matched  materials  systems.  This  capability,  however, 
can  be  increased  substantially  by  including  lattice-matched 
heterostructures  between  chakopyrite  D-IV-VIj  compounds 
and  alloys  [1]  and  the  above  sphalerite  materials.  In  this  letter 
we  present  experimental  evidence  on  and  Hiarw  the  use  of 
p-ZnSnPj  on  n-GaAs  hetero) unctions  to  emit  and  detect 
radiation  below  the  bandgaps  of  the  constituent  senriconduc- 
tors. 

The  band  diagram  for  ZaSaP^GaAs  in  forward  bias  is 
shown  ia  Fig.  1.  This  heterojunctiou  has  an  energy-band 
lineup,  which,  as  discussed  previously  (2],  may  provide  a 
source  of  below  bandgap  emission.  The  mechanism,  indi¬ 
cated  in  Fig.  1,  involves  the  tunneling  of  electrons  across  the 
interfacial  region  either  to  or  from  virtual  states  in  the  gap. 
This  enables  long-wavelength  radiative  recombination  of  the 
confined  electrons  and  holes  in  adjacent  wells  on  the  two 
sides  of  the  interface.  Although  the  band  diagram  in  Fig.  1  ia 
for  a  forward-biased  heterojuaction.  the  equilibrium  configu¬ 
ration  is  very  similar  only  with  deeper  wells  and  wider 
space-charge  regions.  Therefore,  long-wavelength  absorp¬ 
tion  also  appears  possible  involving  the  inveree  of  tha  process 
just  described  (3). 

Photovoltaic,  photoconductivity.  Hail,  I-V,  and  CV  mee- 

pnranwaara  to  coaatruct  an  inmrfhcial  band  diagram.  This 
construction  yielded  a  value  for  the  interfacial  energy  E,.  The 
photovoltaic  and  emiaaioa  experiments  revealed  a  below 
bandgap  peak  which  is  compared  with  the  theoretical  value  of 
E,.  Other  sigaificaat  aspects  of  the  data  aa  they  pertain  to  the 
model  are  diacntasd 
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0.  DEVICE  FABRICATION 

The  heterostructures  used  in  this  investigation  consisted  of 
liquid-phase  epitaxial  (LPE)  layers  of  ZnSnP2  grown  on 
{211}  Ga-GaAs  with  a  lattice  mismatch  of  less  than  0.44 
percent.  The  details  of  the  growth  process  were  reported 
earlier  [4].  The  {211}  orientation  of  the  GaAs  substrate  was 
chosen  since  it  has  been  shown  to  suppress  antiphase  domain 
boundaries  associated  with  different  c-axis  orientations  in  the 
cfaalcopyrite  ZnSnP2  layer  [3}.  The  devices  were  defined 
with  standard  photolithographic  techniques  and  mesa-etched 
with  Br-methanol  to  an  area  of  2  x  10' ’  cm2. 

The  GaAa  substrates  used  were  n*(»  10“  cm'5),  and  the 
ZnSnP2  was  grown  p^fwlO1*  cm~J),  according  to  Hall 
measurements  on  ZnSnP2  layers  grown  on  semi-insulating 
(SI)  GaAs.  There  was  some  evidence  of  Zn  diffusion  into  the 
GaAs  substrates.  Delineation  etches  of  ZnSnP2  layers  grown 
on  SI  substrates  revealed  both  a  growth  interface  and  a 
diffusion  from  of  shoot  1  pm  into  the  substrate.  When  the 
ZnSaP2  growth  was  performed  on  n'-GaAs,  however,  the 
diffosioa  of  Za  was  not  large  enough  to  convert  the  GaAs  to 
p-type,  aa  evidenced  by  a  single  line  visible  at  the  interface. 
In  a  similar  manner,  electron  microprobe  analysis  indicated 
omdiffosioe  of  Ga  (but  not  As)  into  leu  than  1  pm  of  the 
ZaSnP2  with  no  type  conversion. 

ffl.  RESULTS  AND  DISCUSSION 

The  bendgaps  of  the  two  materials  were  obtained  through 
photoconductivity  and  photovoltaic  measurements.  The 
value  for  the  n*-GaAs  bandgap  wax  determined  to  be  1.42 
eV.  The  ZnSnP2  bandgap  of  1.43  eV  was  determined  from 
differential  analysis  of  photoconductivity  data  for  layers 
grown  on  SI  GaAs.  This  method  wu  used  since  a  clear  peak 
for  ZnSaP2  wu  not  apparent  in  the  photovoltaic  spectra  on 
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n*-GaAs.  The  value  of  1.45  eV  for  ZnSuP2  is  consistent  with 
previous  measurements  [6]. 

CV  measurements  were  performed  to  obtain  a  value  for 
the  built-in  potential  K0.  This  method  uses  a  linear  extrapola¬ 
tion  of  C~ 2  versus  V  to  zero.  The  largest  value  of  Kg  obtained 
in  this  manner  was  1.73  V.  We  used  this  value  in  the  analysis 
that  follows  since  the  measured  K0  should  decrease  when 
interface  charge  (of  either  sign)  is  present.  The  true  V0  can 
only  be  obtained  by  correcting  for  the  interface  charge. 

To  compare  with  this  experimental  value,  a  theoretical 
miww  of  the  built-in  potential  can  be  determined  from  the 
ionization  energies,  the  energy  gaps,  and  the  relative 
poaitioos  of  the  electrochemical  potential  in  the  two  materi¬ 
als.  Using  a  dielectric  model  [7],  we  calculated  an  ionization 
energy  of  5.91  eV  for  ZnSnP2  and  used  an  experimental 
ionizatioa  energy  of  5.50  eV  for  GaAs  [8].  In  determining 
the  densities  of  states,  we  took  m,*  ■  0.0665  mg  for  GaAs 
[9]  and  m„*  »  0.39  mo  for  ZnSnP2  [6].  The  carrier 
concentrations  were  n  ■  8  x  10IT  cm'1  for  the  GaAs 
substrate  and p  »  5  x  10"  cm'1  for  the  ZnSnP2.  Using  the 
band  parameters  from  these  numbers  yields  a  built-in 
potential  of  1.74  V.  The  close  agreement  between  this  and 
the  experimental  value  of  1.75  V,  however,  should  not  be 
taken  too  seriously.  The  above  parameters  also  allow  us  to 
calculate  a  value  for  the  interfocial  energy  £/,  which 
designates  the  energy  difference  between  the  bottom  of  the 
electron  well  and  the  top  of  the  hole  well.  The  above  numbers 
yield  a  value  of  £,  -  1.27  eV. 

To  compare  the  above  calculated  value  of  £,  with  our 
experimental  value,  consider  the  results  shown  in  Fig.  2.  The 
higher  solid  curve  is  die  photovoltaic  response  (or  absorp¬ 
tion)  of  the  heterpjunction.  The  lower  dashed  curve  is  the 
emission  spectrum  of  the  heterojunction  in  forward  bias  at  a 
constant  current  density  of  10  A-cm'2.  The  higher  energy 
emission  peak  has  its  maximum  at  the  energy  of  the  half 
power  point  of  the  higher  energy  absorption  peak.  From  their 
energies,  these  peaks  appear  to  be  the  emission  and  absorp¬ 
tion  peaks  for  GaAs.  For  a  band-to-band  transition  one 
expects  this  kind  of  relation  between  the  absorption  and 
emission  peaks. 

The  other  set  of  peaks  at  lower  energy  have  their  maxima 
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at  approximately  the  same  value.  This  would  be  expected  for 
transitions  between  the  wells  shown  in  Fig.  1.  Usually  one 
extracts  the  absorption  edge  from  an  absorption  peak  by 
using  the  half  power  point  energy.  For  the  absorption  and 
emission  peaks  associated  with  the  interfocial  gap,  one  might 
extrapolate  the  steep  part  of  the  tail  down  to  the  zero  line 
since  neither  significant  absorption  nor  emission  is  expected 
below  £/  using  our  model.  The  emission  peak  extrapolates  to 
approximately  1.31  eV  and  absorption  peak  to  about  1.30 
eV.  These  values  are  in  reasonable  agreement  with  our 
theoretically  predicted  value  of  1.27  eV. 

An  obvious  concern  in  analyzing  the  data  is  the  possibility 
that  the  low  energy  peak  may  not  be  due  to  the  mechanism  we 
have  described  but  rather  to  transitions  involving  impurity 
levels  in  the  gap.  The  energy  coincidence  of  the  emission  and 
absorption  peaks  contradicts  an  impurity  level  interpretation, 
since  a  shift  in  these  energies  when  impurities  are  involved  is 
usually  observed  due  to  the  Franck-Condon  principle  [10J. 

Also  of  interest  in  both  spectra  is  the  absence  of  an 
observable  ZnSnP2  peak.  In  the  emission  mode,  minority 
carriers  are  injected  into  both  sides  of  the  heterpjunction  and 
participate  in  the  various  recombination  processes.  The 
supply  of  recombining  species  is  much  more  plentiful  on  the 
GaAs  side  due  to  the  much  larger  majority  carrier  concentra¬ 
tion  of  the  ZnSnP2.  Therefore,  it  is  reasonable  that  the  GaAs 
emission  may  be  the  only  band-to-band  transition  observable. 
In  the  absorption  mode,  the  space-charge  region  on  the  GaAs 
side  is  much  larger  than  that  on  the  ZnSnP2  side.  Since  only 
this  portion  of  the  device  contributes  to  the  absorption  signal, 
again  it  is  reasonable  to  only  see  the  GaAs  peak. 

Finally,  in  Fig.  3  we  see  the  emission  spectra  shift  to 
longer  wavelengths  with  increasing  current  through  the 
device.  The  spectra  shown  are  for  10-,  20-.  30-,  and 
40-A  ■  cm ~ 2  constant  current  density.  This  shift  is  most  likely 
due  to  device  heating.  Double-crystal  X-ray  diffraction 
measurements  indicate  no  difference  in  thermal  expansion 
between  ZnSnP,  and  GaAs  from  room  temperature  to  about 
150*C.  Since  the  lattice  constants  on  both  sides  of  the 
interface  expand  together  at  the  same  rate,  we  expect  to  see 
the  two  emission  peaks  move  at  about  the  same  rate  with 
respect  to  current,  as  in  Fig.  3. 

IV.  CONCLUSIONS 

The  results  on  below  bandgap  emission  in  ZnSnP2/GaAs 
heterojunctions  tend  to  confirm  previous  expectations  [2]  that 
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heterostnictures  with  this  type  of  ‘‘staggered”  band  lineup 
have  potential  application  as  photon  emitters.  The  absorption 
results  are  also  consistent  with  the  below  bandgap  absorption 
recently  observed  [3]  in  Ini  .jrGa,P/GaAs  heterojunctions 
with  similar  band  lineups.  In  addition,  our  results  show  that 
chalcopyrite-on-sphalerite  heterostructures  with  reasonably 
good  interfacial  properties  can  be  attained. 
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Advances  in  technology  have  made  possible  the  fabrication  of  rapidly  varying  heterostructures 
which  hold  the  promise  of  important  applications.  We  develop  a  set  of  approximate  treatments  of 
electron  states  in  a  variety  of  layered  heterostructures.  The  approximations  are  all  based  on  the  con¬ 
cept  of  one-band  generalized  Wannier  functions.  Following  a  discussion  of  the  validity  of  this  rep¬ 
resentation,  we  apply  it  to  an  evaluation  of  the  bound  states  in  a  narrow  quantum  well  in  GaAs, 
which  clearly  demonstrates  the  mixing  of  main  and  satellite  valley  states  as  well  as  the  contribution 
of  evanescent  states,  and  of  the  states  of  a  superlattice  in  a  model  structure  of  up  to  20  quantum 
wells.  As  a  final  example  we  discuss  the  application  of  generalized  Wannier  functions  to  the  match¬ 
ing  of  electronic  states  at  a  heterojunction  between  two  model  band  structures  with  different  effec¬ 
tive  masses,  and  compare  the  formalism  with  alternative  approaches  to  this  problem. 


With  the  development  of  special  growth  techniques 
such  as  molecular  beam  epitaxy  (MBE)  and  metal  organic 
chemical  vapor  deposition  (MOCVD),  tailor-made  sub¬ 
micron  semiconductor  heterostructures  can  now  be 
designed.  An  entirely  new  variety  of  structures,  including 
quantum  wells,  superlattices,  and  modulation-doped 
structures,  can  be  conceived  and  some  have  already  led  to 
successful  device  applications.  Quantum  wells  are  being 
investigated  in  semiconductor  lasers1  and  resonant  tunnel¬ 
ing2  for  submillimeter  radiation.  Potential  applications 
for  superlattices  include  nonlinear  optics  in  relation  to  op¬ 
tical  switching  (bistability),3  the  Bloch  oscillator,4  and  the 
Zener  oscillator3  for  submillimeter  radiation.  An  impor¬ 
tant  application  of  modulation  doping6  is  the  high  elec¬ 
tron  mobility  transistor,  a  low-noise  microwave  amplif¬ 
ier,7  and  a  candidate  for  high-speed  logic.* 

These  semiconductor  heterostructures  are  made  of  suc¬ 
cessive  semiconductor  crystal  layers  grown  on  top  of  each 
other.  Provided  the  lattice  parameters  are  closely 
matched,  the  lattices  of  the  semiconductors  essentially 
cohere  with  a  minimal  perturbation.  Due  to  the  special 
growth  techniques  used,  the  spatial  variation  of  the  semi¬ 
conductor  materials  and/or  of  the  doping  can  be  con¬ 
trolled  so  as  to  occur  in  a  few  lattice  parameters. 

The  physical  properties  of  the  structures  of  these  semi¬ 
conductors  are  the  object  of  extensive  theoretical  and  ex¬ 
perimental  effort  with  respect  to  important  potential  de¬ 
vice  applications,  some  of  which  were  mentioned  above, 
but  most  of  which,  however,  remain  unrealized. 

A  fundamental  property  on  which  all  these  heterostruc¬ 
tures  rely  is  the  spatial  variation  of  the  band  gap  and  the 
conduction  and  valence  band  structures.  In  this  paper,  we 
shall  be  concerned  with  the  development  of  a  consistent 
and  simplified  picture  of  spatially  varying  band  struc¬ 
tures. 

Band  structures  derive  from  the  periodic  nature  of  the 
crystal  potential  and  are  therefore  mathematically  well  de¬ 


fined  only  for  an  infinite  crystal.  However,  in  practice, 
the  concept  of  bands  holds  on  a  microscopic  scale.  This  is 
exemplified  by  the  use  of  band  diagrams  in  classical  de¬ 
vice  theory.  The  total  classical  Hamiltonian  H  can  be 
written  as 

ff(k,r)=£(k)-eK(r)  ,  (1) 

where  £(k)  is  the  conduction  or  valence  band  structure,  e 
the  electron  charge,  and  V(r)  the  electrostatic  potential. 
The  symbol  r  denotes  the  spatial  location  and  k  the  Bloch 
wave  vector  or  quasimomentum.  In  a  band  diagram,  usu¬ 
ally  only  the  bottom  of  the  conduction  band  and  the  top 
of  the  valence  band  are  represented,  which  is  equivalent  to 
specifying  k. 

Equation  (1),  together  with  the  acceleration  theorem,9 
constitute  the  semiclassical  picture  on  the  basis  of  which 
one  can  describe  the  ballistic  motion  of  an  electron  or 
derive  (using  the  Boltzmann  transport  equation)  the  stan¬ 
dard  device  equations.10  In  the  heterostructures  discussed 
above,  the  spatial  dimensions  are  made  sufficiently  small 
for  the  quantization  of  the  electron  states  to  be  signifi¬ 
cant.  One  needs  then  to  revert  to  the  quantum- 
mechanical  form  of  Eq.  (1)  where  the  quasimomentum  is 
replaced  by  the  operator 


and  to  solve  the  Schrodinger  equation, 


Equation  (1),  together  with  (2),  is  often  referred  to  as 
the  Wannier  theorem,  the  derivation  of  which  can  be 
found  in  textbooks  on  semiconductor  theory."  The  Wan¬ 
nier  theorem  relies  on  two  basic  assumptions.  The  first  is 
the  one-band  approximation  which  holds  if  the  perturba¬ 
tion  potential  V(r)  is  not  strong  enough  for  band  mixing 
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to  occur  or  interband  transition  such  as  tunneling  to  take 
place.  Secondly,  the  perturbation  potential  should  vary 
slowly  or  smoothly  with  position  so  as  to  constrain  the 
potential  overlap  matrix  to  a  diagonal  representation. 

These  assumptions  are  adequate  for  the  classical  repre¬ 
sentation  in  which  the  electrons  are  essentially  free  and 
can  be  represented  by  a  wave  packet  ¥(r).  For  the  sub¬ 
micron  heterostructures  of  interest  here,  these  assump¬ 
tions  do  not  hold,  as  the  potential  variations  in 
modulation-doped  or  in  spatially  varying  band-gap  struc¬ 
tures  are  both  large  and  can  occur  in  a  few  lattice  parame¬ 
ters.  In  this  paper  we  shall  propose  a  concept  of  general¬ 
ized  band  structures  for  which  the  Wannier  theorem  can 
be  generalized  and  applied  to  submicron  heterostructures. 

First,  we  shall  rederive  the  Wannier  theorem  in  Sec.  I, 
specializing  it  to  one-dimensional  structures.  Since  the 
submicron  heterostructures  grown  by  MBE  or  MOCVD 
are  layered  structures,  a  one-dimensional  picture  can  be 
implemented,  provided  some  assumptions  are  made.  Such 
one-dimensional  models  are  of  heuristic  interest  as  their 
simplicity  promotes  a  better  insight,  and  they  are  some¬ 
times  used  in  semiquantitative  studies. 

In  Sec.  II  we  discuss  the  extension  of  the  Wannier 
theorem,  using  the  generalized  Wannier  functions,  and 
study  in  Sec.  Ill  two  examples  of  heterostructures.  These 
examples  constitute  both  a  test  and  a  demonstration  of  ap¬ 
plications  of  the  generalized  Wannier  picture.  Finally,  we 
present  in  Sec.  IV  the  application  of  the  generalized  Wan¬ 
nier  picture  to  “true"  heterostructures  and  discuss  the 
generalization  of  the  concept  of  bands  to  spatially  varying 
band  structures. 

I.  THE  ONE-DIMENSIONAL  WANNIER  PICTURE 

In  this  section  we  consider  the  problem  of  a  crystal 
electron  in  an  external  one-dimensional  potential  U.  The 
total  Hamiltonian  of  the  electron  is  given  by 

H=H0  +  U&-t), 

where  H0  is  the  unperturbed  crystal  Hamiltonian  and 
U(x)  is  the  one-dimensional  energy  potential  varying 
along  an  axis  (the  device  axis)  represented  by  the  unit  vec¬ 
tor  3. 

We  assume  the  device  axis  to  be  parallel  to  the  lattice 
vector.  We  only  consider  face-centered -cubic  semicon¬ 
ductor  lattices  for  which  a  lattice  vector  R(n)  can  be  writ¬ 
ten  using  the  orthonormal  basis  of  the  Bravais  lattice  as 
R(n)=4n,  with  n  a  set  of  integers  (ni./tj.Rs),  a*  the  lat¬ 
tice  parameter  of  the  cubic  lattice,  and  A  the  matrix 


2 1;  To]' 

We  shall  denote  as  a  the  modulus  of  the  smallest  lattice 
vector  a  parallel  to  3;  a  is  therefore  the  effective  lattice 
parameter  along  the  direction  3,  and  any  lattice  vector 
parallel  to  3  is  written  R=na,  with  n  an  integer. 

A  reciprocal-lattice  vector  K  is  written  in  the  same 
basis  as  K (/)=2?/,  with  /  a  set  of  integers  (/|,/2,/3)  and  B 
the  matrix 


[-1  1  1 

B  =  ^  1-11, 

0  1  1—1 

so  that  AJB  =2ir.  It  is  easy  to  prove  that  for  cubic  lattices 
there  always  exists  a  reciprocal-lattice  vector  parallel  to  a 
given  direct-lattice  vector.  This  is  equivalent  to  finding 
the  constant  c  such  that  K=cR.  A  possible  solution  is 
c  =8 v/a*2,  which  leads  to 

2  1  1 

/=  1  2  1  n. 

1  1  2 

We  shall  denote  q  as  the  modulus  of  the  smallest 
reciprocal-lattice  vector  q  parallel  to  3  so  that  any  re¬ 
ciprocal  vector  parallel  to  3  can  be  written  K=/q,  with  / 
an  integer.  Since  the  lattice  vector  and  the  reciprocal- 
lattice  vector  satisfy  the  relation 

R(n)-K(/)=2im-/  , 

we  have  that  q  and  a  satisfy 

rq=*aq=2irp  , 

with  p  an  integer.  For  the  <100>  direction 
[R=a(  1,0,0)],  we  have  n=(  — 1,1,1),  and  the  smallest 
reciprocal  vector  is  given  by  l  =  (0,1,1),  so  that  a— a* 
and  p  —2. 

In  the  extended  zone  scheme  the  Bloch-function  solu¬ 
tion  ^(k,r)  of 

H0+(k)=E(kWk)  (3) 

is  a  periodic  function  of  k  along  the  device  direction  d. 
We  can  introduce  a  one-dimensional  quasimomentum  k 
along  the  device  direction  3  defined  by 

with  kL  the  transverse  momentum  perpendicular  to  the 
direction  3.  It  follows  that  a  Bloch  function  along  the 
direction  3  can  be  written  as  <Hkd+klfT)=^(k),  which  is 
a  periodic  function  of  k  with  period  q:  <f>(k  +rtq)=<p(k) 
for  n,  an  integer. 

We  can  now  define,  along  the  direction  3,  a  one¬ 
dimensional  Wannier  function  w(m,r,k:)  as  the  spatial 
Fourier  coefficient  of  the  Bloch  function  <Mk), 

w(m,r,kL  )=(  1/t/q  )  f^^^(kd-hkL,r) 

Xcxp(—ikm2ir/q)dk  , 

(4) 

and  reciprocally  we  have 

^(k3+k1,r)  =  (l/v/q  )  2  u>(m,r,k1)exp(ikm  2n/q)  . 

m 

From  the  orthogonality  property  of  the  Bloch  functions 
and  from  Eq.  (3),  one  easily  verifies  that  these  Wannier 
functions  (WF*s)  form  an  orthogonal  set 

(w{n)  |  w(m))=  f_"W(n,r,kL)*u>(m,r,k1)dx'  =  8„m  , 
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with  x'=d-r,  and  satisfy  the  useful  property 
H0w(m,t,kL)=  2 E(n  -m)w(n,t,kL) , 

m 

with 

E(n)=(l/q)  f'l^Elki+kdtexpi—iknfrr/qMk  . 

The  E(n)  are  recognized  as  the  Fourier  coefficients  of  the 
band  structure  along  the  line  k(k)  parallel  to  the  device 
axis. 

For  a  sufficiently  smooth  and  weak  one-dimensional 
potential  t/(d-r),  no  band  mixing  is  expected  and  the 
transverse  momentum  is  conserved.  One  can  then  expand 
the  electron  state  in  Wannier  functions  of  the  band  con¬ 
sidered  with  an  envelope  /( n,t), 

'V='2lf(n,t)w(n,t,kL).  (5) 

* 

The  weight  function  |  f(n,t)  \ 1  is  the  probability  of  find¬ 
ing  the  electron  at  the  lattice  site  n  at  time  t.  Replacing 
¥  in  the  Schrodinger  equation 

HV=i*4~'i>  (6) 

dr 

by  the  expansion  of  Eq.  (3),  multiplying  Eq.  (6)  by  win)*, 
and  integrating  over  x\  one  obtains  the  Wannier  re¬ 
currence  equation 

2  [&in  —m)+  Uin,m)]fim,t)=Hl-^f(n,t)  ,  (7) 

where  U(n,m)  is  the  matrix  element  given  by 

Uin,m)=(win)  |  U  |  ui(m)> 

=  f  win,r,kL)* U(Z-t)wim,r,kL)dx'  .  (8) 

The  one-band  Wannier  recurrence  equation  derived  holds 
for  sufficiently  small  perturbation  potentials  U.  For  Eq. 
(7)  to  be  solved,  both  the  band-structure  Fourier  coeffi¬ 
cients  and  the  matrix  elements  Uin.m )  are  required. 

As  expressed  by  Eq.  (8),  the  matrix  elements  are 
evaluated  from  the  Wannier  functions.  Techniques  for 
the  evaluation  of  the  Wannier  functions  have  been 

developed  by  Kohn.12  The  evaluation  of  the  WF  for  sil¬ 
icon  has  been  recently  reported  by  Kane  and  Kane;13  how¬ 
ever,  the  WF  derived  is  a  linear  superposition  of  the  WFs 
of  the  four  valence  bands.  The  evaluation  of  realistic 
WFs  is  a  complicated  matter,  which  we  would  like  to 
avoid  in  the  simple  picture  we  intend  to  develop. 

Some  direct  approximations  can  be  derived  from  the 
properties  of  the  WFs.  In  the  case  of  simple  bands,  it  is 
possible  to  select  the  phase  of  the  Bloch  functions  so  that 
the  WFs  are  exponentially  localized.14  From  both  the 
orthogonality  and  the  tight  localization  of  the  WFs,  it 
follows  that  the  matrix  element  of  a  smooth  potential  is 
accurately  represented  by  the  sampling  of  the  potential  at 
the  lattice  sites: 

Uin,m)  =  Ulx'  =  na/p)bHm  .  (9) 

We  shall  discuss  in  Sec.  Ill  the  extent  to  which  these  ap¬ 
proximations  can  be  used  for  sharply  varying  potentials. 


Finally,  an  important  property  of  the  one-dimensional 
WFs  which  we  did  not  address  so  far  is  their  invariance 
under  translation  from  one  lattice  site  to  another.  Follow¬ 
ing  the  approach  of  Wannier,13  we  first  write 

♦<k,r— R)=(l/v^  . 

m 

(10) 

which  using  the  translation  property  of  the  Bloch  func¬ 
tion,  can  be  written 

♦(k,r— R)=(l/V^  >2exp<  —  /k-R)exp(ikma/p) 

HI 

Xurim.r.kJ  . 

For  R=na,  we  have 

k-K=ikZ+kl)-K=kZ-K=kna  , 
and  changing  the  index  of  summation  we  obtain 
<Kk,r— ID^il/y/q  )'2>lcxp(ikma/p)w(m  +np,T,kL)  . 

m 

(in 

Identifying  the  coefficients  in  the  series  (10)  and  (11),  we 
finally  get 

w(m,r—nn,kL)—iv<m  +  np, r.k^  . 

Considering  the  <  100)  direction  for  which  p  =2,  it  fol¬ 
lows  that  there  exist  two  sets  of  WFs,  each  of  which  can 
be  generated  from  two  generic  WFs  by  a  lattice  transla¬ 
tion 

wi2n,t)=*wi0,r— /»*),  u>(l+2n,r)=u>(l,r— ns)  . 

The  number  of  WFs  for  a  given  direction  is  therefore 
given  by  p  for  the  face-centered-cubic  lattice  considered. 
For  our  purpose  it  is  convenient  to  picture  these  WFs  as 
located  at  p  sites  along  the  lattice  vector  a  associated  with 
the  direction  3.  We  can  then  write  these  WFs  as 

win,r,kL)=wix’  — na/p,ri,kl)  =  w(x'  —  na/p)*=w(n)  , 

where  a  is  the  effective  lattice  parameter  for  the  direction 
d.  The  one-dimensional  WF  we  have  introduced  is  in  fact 
a  hybrid  entity,  since  it  corresponds  to  the  Bloch  state  kL 
for  the  transverse  coordinates  tL.  The  properties  of  the 
WFs  of  one-dimensional  lattices  cannot  be  assumed  to 
apply  directly  to  the  one-dimensional  WF  we  have  de¬ 
fined.  Therefore,  problems  concerning  the  reality  of  the 
WF  and  the  simultaneous  convergence  of  the  WFs 
remain  to  be  addressed.  These  studies  are  not  trivial 
matters  for  three-dimensional  crystals,14  in  particular  for 
lattices  without  a  center  of  inversion. 

In  our  simplified  treatment  we  shall  assume  the  picture 
developed  to  be  truly  one-dimensional.  The  reality  and 
exponential  localization  of  the  WFs  in  one-dimensional 
structures  are  discussed  by  Kohn.16  More  recent  work  by 
Zak17  raises  the  one-dimensional  Wannier  picture  to  the 
status  of  a  well-defined  representation  in  quantum 
mechanics  and  introduces  the  canonical  Wannier  function 
with  minimal  position  uncertainty  and  located  in  the  zero 
cell  of  the  Bravais  lattice. 

Zak  also  demonstrates  that  the  position  operator  has  a 
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discrete  spectrum.  This  implies  that  the  applied  poten¬ 
tials  are  transformed  by  the  one-band  approximation  into 
steplike  potentials.  This  is  due  to  the  inability  of  the  one- 
band  approximation  to  deal  with  the  perturbation  of  the 
inner  crystal  potential.  However,  as  we  shall  see  in  the 
next  section,  most  of  the  band  mixing  can  be  formally  ac¬ 
counted  for  with  the  use  of  the  generalized  Wannier  func¬ 
tion.  For  stronger  potentials  it  is  further  necessary  to  ac¬ 
count  for  the  polarization  of  the  semiconductor.11 

II.  THE  GENERALIZED  WANNIER  PICTURE 

In  this  section  we  shall  be  concerned  only  with  abrupt 
modulation-doped  heterostructures  or  models  of  abrupt 

- 1 


heterostructures  for  which  the  band  gap  is  spatially  vary¬ 
ing  while  the  band  structure  is  assumed  to  be  conserved. 

The  one-dimensional  Wannier  picture  developed  in  the 
previous  section  is  capable  of  handling  fast  spatial  varia¬ 
tion  of  the  potential  since  the  nondiagonal  term  of  the 
matrix  potential  V(n,m)  can  be  accounted  for.  However, 
the  formalism  developed  relies  on  the  one-band  approxi¬ 
mation  which  strictly  holds  only  for  small  potential  per¬ 
turbations.  Its  use  for  submicron  heterostructures  seems 
unjustifiable.  Indeed,  the  band-edge  variation  in  these 
structures  is  large  and  band  mixing  is  expected. 

One  should  then  resort  to  a  multiband  expansion,  and 
the  envelope  equation  (7)  now  reads 


where  b  is  the  band  index. 

It  is,  however,  possible  to  avoid  a  multiband  expansion 
if  one  introduces  the  concept  of  generalized  bands.  As 
was  remarked  by  Zak,17  “a  band  is  a  surprisingly  stable 
entity  under  external  perturbations.”  This  might  seem  to 
be  a  paradoxical  statement  coming  from  7-ak,  who  has 
been  very  critical  of  the  use  of  the  one-band  approxima¬ 
tion  in  general,19  and  in  particular  for  the  Stark  effect.30 
In  fact,  the  existence  of  bands  under  external  perturba¬ 
tions  does  not  rely  on  the  one-band  approximations.  A 
band  of  energy  is  essentially  a  quasicontinuum  of  energy 
states.  Such  an  entity  seems  to  survive  numerous  mis¬ 
treatments  such  as  high  impurity  concentrations,  alloying 
with  random  distribution  of  two  atoms  on  a  given  sublat¬ 
tice  in  III-V  compounds,  even  for  abrupt  heterqjunctions. 
However,  surface  states  or  band  tailing  or  even  alloy  and 
roughness  scattering  effects  occur,  and  it  might  be  more 
accurate  to  refer  to  these  bands  as  generalized  bands. 
Such  bands  in  most  cases  can  only  be  represented 
mathematically  as  resulting  from  the  mixing  of  the  bands 
of  the  unperturbed  crystals.  This  is  similar  to  the  expan¬ 
sion  of  the  pseudo-wave-function  in  terms  of  plane  waves 
in  the  calculation  of  band  structures  using  the  pseudopo¬ 
tential  method;  the  bands  evaluated  result  from  the  mix¬ 
ing  of  the  bands  of  the  empty  lattice. 

A  one-band  formalism  can  then  be  developed  if  one 
uses  the  exact  Wannier  functions  instead  of  the  Wannier 
functions  of  the  unperturbed  lattice.  It  remains  to  be 
demonstrated  that  such  Wannier  functions  exist.  The 
Wannier  functions  as  we  introduced  than  [see  Eq.  (4)]  are 
defined  as  Fourier  coefficients  of  the  Bloch  functions 
which  are  periodic  in  k  space.  The  quasimomentum  is 
defined  by  the  translation  operator11  only  for  periodical 
potentials.  At  a  heterojunction,  the  periodicity  is  broken 
and  no  quasimomentum,  Bloch  functions,  or  Wannier 
functions  can  be  defined,  but  one  can  attempt  to  general¬ 
ize  the  entities,  as  the  periodicity  is  essentially  broken  only 
at  the  junction.  The  most  obvious  approach  is  to  general¬ 
ize  the  Wannier  functions  which  are  strongly  localized  on 
each  lattice  site  and  therefore  will  be  strongly  perturbed 
only  at  the  junction  itself. 


The  existence  of  generalized  Wannier  functions 
(GWFs)  has  been  theoretically  established  by  Kohn  and 
Onffroy21  for  one-dimensional  structures.  These  Wannier 
functions  can  be  labeled  by  their  generalized  band  index; 
they  account,  however,  for  the  presence  of  interface  states 
in  the  forbidden  band  gap.22  We  shall  denote  these  gen¬ 
eralized  Wannier  functions  as  a(n,x).  They  span  the 
same  function  space  as  the  eigenfunctions  of  the  per¬ 
turbed  lattice  for  the  same  generalized  band. 

The  index  n  will  be  taken  here  as  the  site  index,  al¬ 
though  depending  on  the  potential  well  there  may  be  more 
than  one  or  even  no  GWFs  associated  with  the  interface 
sites.  The  GWFs  are  localized  about  the  lattice  site  (or 
the  periodic  potential  minima).  Rehr  and  Kohn  have 
shown23  that  they  can  be  exponentially  localized.  An 
essential  property  of  the  GWFs  is  that  they  form  an 
ortkonormal  set 

(a(n)\a(m))=  a(n,x)a(m,x)dx  =  8m  . 

The  eigenfunctions  are  linear  combinations  of  a  (n,x), 
'Kx,E)=  2f(n,E)a(n,x) , 

a 

with  /( n,E)  the  solution  of 

£/(«,£)=  2mn,m)f(m,E)  ,  (12) 

m 

where  H(n,m)  is  the  Hamiltonian  matrix  element  given 
by 

H(n,m)=(a(n)\H \a(m))  . 

Here  H  is  the  total  Hamiltonian  of  the  heterojunction 
given  by 

H=H0  +  U(x)  , 

with  H0  the  unperturbed  crystal  potential  and  U(x)  the 
heterojunction  effective  potential. 

Equation  (12)  is  similar  to  Eq.  (7),  except  that  we  use 
the  GWF  instead  of  the  unperturbed  Wannier  functions. 
The  knowledge  of  the  GWF  weight  /( n,E)  is  sufficient  in 
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most  applications  and  Eq.  (12)  can  be  solved  if  the  matrix 
dement  Hln,m)  is  known. 

The  derivation  of  the  matrix  dement  is  the  critical  part 
of  this  method.  One  approach  consists  of  first  deriving 
the  GWFs  at  each  site.  An  important  property  of  the 
GWFs  demonstrated  by  Kohn  and  Onffroy  states  that 
for  lattice  sites  far  from  the  interface  the  GWFs  ap¬ 
proach  the  unperturbed  WF  exponentially  in  it, 

a(n,x)stw(x  — na )  for  | n  |  »0  , 

assuming  the  interface  at  n  =0.  Only  a  few  GWFs  are 
found  to  be  perturbed  or  to  differ  from  the  WF,  as  was  re¬ 
ported  by  Gay  and  Smith.22  It  follows  that  a  few  lattice 
parameters  away  from  the  interface  the  Wannier  picture 
devdoped  in  Sec.  I  holds  and  the  only  remaining  unknown 
terms  are  the  interface  Hamiltonian  matrix  terms. 

A  variational  method  for  the  derivation  of  the  GWF 
has  been  proposed  by  Kohn  and  Onffroy  for  one¬ 
dimensional  potentials  and  was  tested  by  Gay  and 
Smith.22  A  three-dimensional  treatment  is  required  for 
the  heterostructures  of  interest  here.  The  derivation  of 
the  Wannier  functions  in  three  dimensions  being  quite  in¬ 
volved,11  the  derivation  of  the  more  complicated  GWF 
will  not  be  considered  as  a  practical  approach.  For  our 
present  purposes,  we  shall  use  for  the  Hamiltonian  matrix 
simple  approximations  derived  from  the  general  proper¬ 
ties  of  the  Wannier  functions  discussed  in  Secs.  I  and  II. 

The  essential  result  brought  about  by  the  GWF  formal¬ 
ism  is  that  as  long  as  a  generalized  band  of  energy  can  be 
assumed,  it  is  possible  to  lump  all  the  band-mixing  effects 
into  the  interface  Hamiltonian  matrix  dements,  whereas 
the  wave  function  fin)  remains  given  by  an  equivalent 
Hamiltonian  recurrence  equation. 

III.  TWO  SIMPLE  EXAMPLES 

We  shall  now  consider  two  simple  applications  of  the 
Wannier  picture.  The  first  consists  of  the  evaluation  of 
the  resonant  states  in  a  square  quantum  well  0.3  eV  deep 
using  GaAs  band  structure.  One  expects  to  observe  the 
interaction  of  the  upper-valley  and  central-valley  states  of 
GaAs  for  such  a  deep  well.24 

These  band-structure  data,  kindly  furnished  by  Profes¬ 
sor  Karl  Hess  of  the  University  of  Illinois,  were  given  at  a 
set  of  sampled  points  in  the  Brillouin  zone.  Using  the 


FIG.  1.  GaAs  band  structure  along  A  (dots:  original  sampled 
data;  solid  curve:  reconstructed  band  structure  with  effective 
mass  fitted  at  T  to  0.068m). 
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FIG.  2.  Fourier  coefficients  of  the  T  to  X  band  structure  of 
Fig.  1. 


sampling  theorem,  we  reconstituted  the  band  structure. 
The  band  structure  was  in  addition  fitted  at  T  to 
0.068  m0.  This  consequently  slightly  altered  the  smooth¬ 
ness  of  the  upper  valley.  The  resulting  band  structure  of 
GaAs  along  A  is  shown  in  Fig.  1.  The  amplitudes  of  the 
Fourier  coefficients  of  the  GaAs  band  structure  along  A 
are  plotted  in  Fig.  2.  We  present  in  Fig.  3  the  locus  of  the 
amplitude  of  the  roots  exp \jki(E)a]  that  solve  the  equa¬ 
tion  E(kt)—E,  where  Elk)  is  the  band  structure  along  A. 
One  obtains  20  roots.  Roots  1  and  2  correspond  to  Bloch 
wave  functions  of  the  central  valley  at  T.  Roots  3  and  4 
correspond  to  the  Bloch  waves  of  the  upper  valley  cen¬ 
tered  on  X.  The  unexpected  behavior  at  X  of  roots  3, 6,  7, 
and  8  is  due  to  the  fact  that  the  upper-valley  minimum 
does  not  occur  exactly  at  X.  The  shift  of  the  X  minimum 


FIG.  3.  Locus  of  the  roots  of  GaAs  band  structure  of  Fig.  1 
as  a  function  of  the  energy. 
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is  an  artifact  of  the  spectrum  truncation.  Since  the  band 
structure  is  slightly  distorted,  we  do  not  expect  these  roots 
to  contribute  much.  There  are  finally  12  remaining  roots 
with  a  very  high  damping  rate.  These  roots  contribute 
only  within  the  matching  area.  One  observes  in  Fig.  3  the 
evanescent  Bloch  waves  being  transformed  into  propa¬ 
gating  Bloch  waves  at  T  (0  eV)  and  X  (0.3  eV). 

The  Wannier  picture  calls  for  the  matrix  element  of  the 
quantum  well  potential.  For  trial  purposes  we  evaluated 
numerically  (with  the  Gauss- Legendre  integration  tech¬ 
nique)  its  matrix  element  using  the  WFs,  the  weakly  lo¬ 
calized  WFs  of  the  empty  lattice 
1 1  n 


w(x  —ma/p)= 


sin  Upx/a  —mhr] 
[(px/a  —m)  v\ 


It  is  most  enlightening  to  consider  the  diagonal  terms 
plotted  in  Fig.  4  (with  an  arbitrary  energy  scale).  One  ob¬ 
serves  that  these  diagonal  terms  approximately  follow  the 
assumed  square  potential,  even  for  such  a  narrow  well. 
The  energy  difference  (a  few  percent)  is  distributed  in  the 
nondiagonal  terms  and  contributes  somewhat  to  the 
matching.  Using  the  complete  matrix  element,  we  found 
that  sufficiently  accurate  results  could  be  generated  by 
directly  reducing  the  matrix  dement  to  diagonal  terms 
corresponding  to  the  sampling  of  the  assumed  potential  at 
the  lattice  site  [see  Eq.  (9)]. 

We  have  devdoped  an  algorithm  to  solve  the  eigenvalue 
problem  and  have  evaluated  the  locus  of  the  resonant  en¬ 
ergy  levd  of  the  square  quantum  wdl  of  depth  0.5  eV,  as 
a  function  of  the  well  width  L.  The  result  is  plotted  in 
Fig.  5  (open  circles)  and  compared  with  the  energy-levd 
locus  (solid  curve)  of  the  effective-mass  approximation.  It 
should  be  pointed  out  that  in  our  picture  the  wdl  width 
can  only  be  varied  by  increments  of  a  half-lattice  parame¬ 
ter  for  the  100  direction. 

We  labd  the  first,  resonant  state  1.  For  a  wide  poten¬ 
tial  well,  both  the  WF  and  the  continuum  methods  yidd 
the  same  results,  since  the  resonant  levd  lies  at  low  energy 
where  the  effective  mass  is  quite  accurate.  For  a  narrower 
wdl,  one  observes  a  divergence  of  the  two  energy  plots  as 
the  effective-mass  approximation  becomes  inadequate  at 
higher  energy. 

For  a  wdl  smaller  than  five  lattice  parameters,  the  first 
resonant  state  occurs  above  0.3  eV.  This  resonant  state  re¬ 
sults,  then,  from  a  superposition  of  the  central-valley 
Bloch  wave  and  the  upper-valley  Bloch  wave.  Indeed,  one 
observes  the  interference  of  the  energy  locus  (labded  a)  of 


FIG.  5.  Energy  levels  of  a  GaAs  quantum  well,  as  a  function 
of  the  well  width.  The  well  depth  is  0.5  eV. 

the  first  resonant  levd  of  the  upper  valley  with  the  energy 
locus  1  of  the  central  valley.  This  destroys  the  smooth 
variation  of  these  energy  lod  as  they  interact.  The  same 
phenomenon  occurs  also  as  the  second  and  third 
resonant-levd  lod  (denoted,  respectively,  2  and  3)  cross 
the  X  relative  minimum  (0.3  eV).  A  more  thorough  dis¬ 
cussion  of  this  interference  effect  has  been  recently  report¬ 
ed  by  Chang  and  Ting.24 

The  phenomenon  just  described  constitutes  an  example 
of  what  may  be  called  a  full  band-structure  effect  which 
can  be  handled  by  the  generalized  Wannier  picture.  The 
algorithm  devdoped  can  be  used  to  study  any  type  of 
variation  of  the  band  gap.  As  a  second  example,  we  con¬ 
sider  the  simple  case  of  a  superlattice  made  of  square 
wells  such  as  in  the  Kronig-Penney  problem.25  We  as¬ 
sume  a  simple  tight-binding  band  structure 
E(k)=A  (1—  coakfl)  with  the  effective  mass  at  k—0 
selected  to  be  0.12m0.  The  wdl  depth  is  0.2  eV,  the  width 
and  separation  of  the  wells  is  chosen  to  be  10a  so  that  the 
superlattice  parameter  is  20a.  This  leads  to  two  resonant 
levels  in  an  individual  wdl.  In  Fig.  6  we  present  the  result 
for,  successivdy,  1,  2,  and  20  wells.  The  discrete  band 
structure  of  the  20-well  superlattice  is  plotted  in  Fig.  7 
(open  circles).  This  is  to  be  compared  to  the  band  struc¬ 
ture  (solid  curve)  of  the  Kronig-Penney  model  in  the 
effective-mass  approximation. 


FIG.  4.  Diagonal  matrix  element  of  a  quantum  well  Uix), 
nine  lattice  parameters  wide. 
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FIO.  7.  Superlattice  band  structure  (open  circles:  tight* 
binding  superlattice  of  20  quantum  wells;  solid  curve:  Kronig- 
Penney  model). 


The  low-lying  bands  (band  1)  have  a  similar  shape  but 
are  separated  by  a  shift  of  1  meV.  This  small  shift 
originates  from  the  use  of  the  approximate  diagonal  ma¬ 
trix  dement  for  the  superlattice  potential. 

The  upper  bands  (band  2),  however,  differ  appreciably. 
Part  of  the  difference  is  due  to  the  deviation  from  para- 
bolicity  of  the  cosine  band  and  results  in  a  much  higher 
conduction-band  maximum  for  the  effective-mass  approx¬ 
imation.  The  remaining  difference  is  essentially  due  to 
the  finite  character  of  the  superlattice  composed  of  only 
20  wells. 

More  sophisticated  structures  can  be  studied  with  the 
algorithm  developed,  but  this  was  not  our  purpose,  as  we 
intended  to  test  the  generalized  Wannier  picture  using 
simple  examples  for  which  closed-form  solutions  are 
available. 


IV.  SPATIALLY  VARYING  BAND  STRUCTURES 

The  Wannier  picture  devdoped  above  and  its  generali¬ 
zation  so  far  do  not  apply  to  band  structures  varying  with 
space.  Let  us  point  out  right  away  that  the  ides  of  spa¬ 
tially  varying  band  structures  is  an  ill-defined  concept, 
which  can  be  made  more  rigorous  within  the  framework 
of  the  generalized  Wannier  picture. 

A  one-dimensional  approach  will  break  down  for  the 
proper  description  of  a  spatially  varying  band  structure, 
since  the  variation  in  the  band  structure  occurs  in  the 
three-dimensional  k  space  and  therefore  requires  a  three- 
dimensional  matching  of  the  Wannier  envelope. 

We  shall,  however,  for  simplicity  partially  retain  the 
one-dimensional  picture  by  assuming  the  spatially  varying 
band  structure  to  be  given  by 


#k\ 

E(k,x)=*E(k,x)+- — j-  , 

2m* 

where  m*  is  a  transverse  effective  mass  which  is  not 
varying  in  space.  Such  a  band  structure  does  not  corre¬ 
spond  to  any  realizable  interface  but  constitutes  a  con¬ 
venient  model.  In  this  modd  we  use  the  one-dimensional 
GWF  formalism  to  solve  for  the  eigenstates  of  this  struc¬ 
ture.  A  set  of  orthogonal  GWFs  noted  a  (n,x)  can  be  de¬ 
fined  as  in  Sec.  n.  These  GWFs  can  be  exponentially  lo¬ 
calized  at  each  lattice  site  n.  For  the  case  of  an  abrupt 
heterpjunction,  the  GWFs  a(n,x)  approach  exponentially 
with  n  the  WFs  of  semiconductors  1  and  2,  for  lattice 
sites  far  from  the  interface  (at  n  =0).  A  few  lattice  sites 
away  from  the  interface,  inside  semiconductor  1,  a  GWF 
closdy  resembles  the  WF  of  the  unperturbed  semicon¬ 
ductor  1;  the  orthogonality  with  the  GWFs  of  semicon¬ 
ductor  2  across  the  interface  essentially  affects  only  the 
exponentially  decaying  tail  of  this  GWF.  Finally,  one  is 
again  led  to  the  generalized  one-band  formalism  of  Sec.  II, 
with  the  GWF  envdope  given  by  the  one-dimensional 
Wannier  recurrence  equation  (7). 

The  derivation  of  the  Hamiltonian  matrix  dements 
H(n,m)  is  again  critical  to  the  analysis.  We  favor  a 
direct  evaluation  of  these  dements,  which  remains  to  be 
devised.  In  its  absence,  alternative  approaches  can  be  con¬ 
caved,  and  we  shall  demonstrate  such  an  approach  for  a 
simple  system. 

First  consider  the  trial  Hamiltonian 

H(n,m)=E[n  —m,n)  +  U(n,m) , 

where  E(h  — m,n )  is  the  band-structure  Fourier  coeffi¬ 
cient  which  is  assumed  to  vary  in  space  with  the  lattice 
site  n.  Such  a  Hamiltonian,  which  might  be  proposed  for 
a  smoothly  spatially-varying  band  structure,  is  not  correct 
as  the  Hamiltonian  is  no  longer  Hennitian: 

H(n,m)^zH*(m,n)  . 

This  arises  as  a  consequence  of  using  the  intrinsically  ill- 
defined  concept  of  spatially  dependent  band  structures. 
The  GWF  formalism  is  capable  of  accounting  for  the 
breakdown  of  the  band  structure.  The  Hermiticity  is 
indeed  compatible  with 

E(m  — Rt«)#£(n  —  m,n)  ; 

it  follows  that  the  related  band  structure  at  the  lattice  site 
n, 

E(k,n)=2lE{l,n)el*k , 

i 

is  now  imaginary  and  therefore,  as  expected,  has  lost  its 
usual  physical  meaning. 

For  an  abrupt  heterojunction,  the  use  of  GWF  formal¬ 
ism  is  reduced  to  the  derivation  of  matching  rules. 
Indeed,  away  from  the  interface  the  Wannier  picture  of 
Sec.  I  holds  and  the  only  unknown  terms  are  the  interface 
Hamiltonian  element  H(n,m). 

Consider  the  tight-binding  band  structures  given  by 

E\(k,x)=*A  —A  cos (ka),  Ei(k,x)=B—B<xx(ka) , 
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where  the  amplitudes  A  and  B  are,  respectively,  related  to 
their  effective  mesa  by 

.  d2  ,  _  # 

A  - - -  and  B=* - r  . 

m^a 

For  simplicity  we  evaluate  the  matrix  dement  H(n,m) 
using  the  unperturbed  WFs  on  both  sides  of  the  junction. 
In  the  tight-binding  approximation,  the  only  remaining 
unknown  is  then  the  matrix  element  £f(i,i  +  1)»— C/2, 
assuming  the  interface  is  located  somewhere  between  the 
lattice  site  i  and  i  +  l.  The  resulting  Hamiltonian  is 
shown  around  the  interface  in  Fig.  8  for  l/(«,m)=0. 

The  evaluation  of  C  depends  now  on  the  particular 
matching  theory  upon  which  we  choose  to  rely.  There  has 
been  much  published  recently  on  the  matching  of  wave 
functions  across  a  heterojunction.  The  GWF  picture  pro¬ 
poses  a  general  method  which  enables  us  to  evaluate  the 
wave  function  in  heterostructures  but  calls  for  the 
knowledge  of  H(n,m).  In  its  absence  we  intend  first  to 
relate  this  method  to  other  matching  theories. 

Popular  techniques  are  the  effective-mass  matching  (see 
Ref.  26  and  the  discussion  by  Kroemer  and  Qi-Gao 
Zhu,27,2*  matching  developed  from  the  (k-p)  Kane  model 
by  Bastard,2*  and  White,  Margues,  and  Sham,30,31  and 
from  tight-binding  arguments  by  White,  Margues,  and 
Sham30  and  Ando  and  Mori.32  A  review  of  the  different 
methods  is  briefly  given  in  White,  Margues,  and  Sham.30 
We  intend  to  give  a  more  complete  discussion  of  the 
matching  problem  in  a  forthcoming  paper. 

It  follows  from  the  work  of  Zhu  and  Kroemer2*  that 
the  effective-mass  matching  is  related  for  small  energies 
to  the  geometrical  avenge  for  the  “ideal”  heterostructure 
(type  L  Ref.  28)  considered  here, 

C=(AB)'n  . 

From  Bastard’s  work,29  we  have  C  given  by  the  average 
CMA+B)/ 2  . 

Both  expressions  actually  lead  to  the  same  results  when 
the  variation  of  effective  mass  is  small  across  the  hetero¬ 
junction.  In  Fig.  9  we  compare  the  effective- mas*  match¬ 
ing  with  the  geometrical  average  matching  for  a  quantum 
well.  This  quantum  well  is  made  of  a  layer  of 
Gao.47Ino.5jA*  sandwiched  between  Alo.4jIno.52Aa.  The 
effective  masses  are,  respectively,  taken  as  0.041  and 
0.073  m0  for  GalnAs  and  AlInA*.  These  effective  masses 
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FIG.  8.  Tight-binding  Hamiltonian. 


FIG.  9.  Energy  locos  as  a  function  of  well  width  for 
Gaa47liVi  jjA*/ Alo. wina  5A1  quantum  well  (open  circles:  results 
using  Wannier  algorithm  with  the  geometrical  average  match¬ 
ing;  solid  curve:  results  using  the  effective  mass  matching).  The 
well  depth  is  0.3  eV. 


should  be  corrected  for  narrow  wells,  as  the  nonpnraboii- 
city  is  no  longer  negligible  (see  Welch,  Wicks,  and  East¬ 
man33).  The  discrete-energy  locus  is  seen  to  agree  with 
the  continuous-energy  locus  baaed  on  the  effective-mass 
formalism. 

This  result  eras  established  within  the  tight-binding  ap¬ 
proximation.  A  more  refined  scheme  will  extend  the  con¬ 
nection  rules  to  include  higher-order  overlaps. 


V.  PERSPECTIVE 

We  have  introduced  a  generalized  one-dimensional 
Wannier  picture.  This  formalism  enables  us  to  handle 
sharp  variations  in  potential  as  long  as  the  potential 
strength  is  sufficiently  moderate  for  the  concept  of  a  gen¬ 
eralized  band  to  bold.  The  associated  band-mixing  effects 
are  then  lumped  into  a  generalized  one-band  Hamiltonian 
matrix.  In  the  absence  of  a  knowledge  of  the  GWF,  one 
can  often  rely  in  practice  on  the  direct  sampling  of  the  po¬ 
tential  at  each  lattice  site. 

The  generalized  one-band  picture  applies  directly  to 
modulation-doped  structures  and  can  be  used  as  a  model 
of  heterostructures  by  trading  the  band-gap  variation  as 
an  effective  potential.  Another  application  not  covered 
here  and  for  which  this  model  was  in  fact  developed,  con¬ 
cerned  the  study  of  devices  such  as  the  Zener  superlattice 
oscillator,  for  which  the  periodicity  in  k  space  of  the  sub¬ 
band  together  with  its  spatial  variation  are  fundamental  to 
the  device  operation.3  Finally,  practical  derivation  of  the 
energy  levels  and  eigenfunctions  can  be  performed  for  ar¬ 
bitrary  potentials  or  heterostructure  configurations  using 
the  same  algorithm. 

The  one-dimensional  generalized  Wannier  formalism 
can  be  used  to  model  true  heterostructures,  although  a 
rigorous  treatment  requires  a  three-dimensional  GWF 
analysis.  In  the  one-dimensional  picture,  the  matching  of 
the  wave  function  across  the  interface  is  determined  by 
the  interface  Hamiltonian  matrix  elements. 
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The  derivation  of  Wannier  functions  is  in  general  a 
complex  problem  due  to  their  orthogonality  and  the  fact 
that  they  are  not  eigenfunctions  of  the  Hamiltonian.  The 
best  procedure  would  be  to  evaluate  the  generalized 
Fourier  coefficients  of  the  generalized  band  structure 
directly,  just  as  the  band  structure  is  evaluated  directly  in 
pseudopotential  calculations  using  variational  techniques. 

An  alternative  approach  is  to  derive  the  interface  Ham¬ 
iltonian  matrix  directly  from  the  enforcement  of  physical 
constraints.  Such  an  approach  was  undertaken  by  Zhu 
and  Kroemer2*  and  Ando  and  Mori32  for  maximized 
matching  (ideal  heterojunction)  using  current  continuity. 
We  shall  report,  in  a  forthcoming  paper,  on  the  extension 
of  these  matching  theories  to  a  higher  energy  range  (non¬ 
effective-mass  case)  and  for  the  inclusion  of  higher-order 
overlap  (non-tight-binding  cases)  for  both  type  I  and  type 
II  semiconductors  (see  Ref.  28). 

The  concept  of  generalized  bands  described  here  is 
largely  intuitive  and  has  been  used  so  far  as  such.  The  in¬ 
troduction  of  the  generalized  Wannier  functions  brings 
support  to  the  concept.  One  interesting  feature  of  the 
generalized  Wannier  picture  is  that  it  accounts  consistent¬ 
ly  for  the  concept  of  spatially  varying  band  structure. 
The  generalized  Wannier  picture  appears  then  as  an  inter¬ 
mediate  method  situated  between  the  direct  study  at  the 
lattice  level  [e.g.,  linear  combination  of  atomic  orbitals 
(LCAO)  techniques]  and  the  techniques  based  on  the 
effective-mass  approximation. 


We  believe  that  useful  insights  can  be  drawn  from  this 
ability  of  the  GWF  formalism  to  handle  the  concept  of 
band  structure,  through  the  use  of  the  band-structure 
Fourier  coefficients  together  with  their  spatial  dependence 
on  the  lattice  sites.  The  band  structure  appears  as  a  phys¬ 
ical  entity  established  in  a  few  lattice  parameters.  We 
found  that  five  to  ten  lattice  sites  were  sufficient  for  a 
close  representation  of  a  band  structure  along  <100). 
This  is  in  qualitative  agreement  with  the  Wannier  func¬ 
tions  evaluated  by  Kane  and  Kane.13  Energy  bands  calcu¬ 
lated  from  these  Wannier  functions  were  accurate  to  0.1 
eV  when  the  Wannier-function  overlaps  were  evaluated  to 
twelfth  neighbors  (439  individual  bonds).  The  spatial  es¬ 
tablishment  of  a  small  effective  mass  requires  a  larger 
number  of  sites  since  a  larger  number  of  Wannier- 
function  overlaps  or  band-structure  Fourier  coefficients  is 
required  to  fit  a  small  effective  mass  while  accurately 
reproducing  the  remainder  of  the  band  structure. 
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10.  LONG  RANGE  ORDERING 

For  a  semiconductor  alloy,  AxBi_xc*  regular  solution 
theory  indicates  that  the  excess  enthalpy  or  heat  of  mixing 
is  given  by 

a  Hm  =  x(l-x) n, 

where  a  is  an  interaction  parameter  which  contains  the  physics 
of  the  mixing  process.  For  positive  values  of  n  we  expect 
to  have  cation  clustering  above  and  phase  separation  below 
a  critical  temperature. 


Tc  =  k41n2  * 

For  negative  values  of  Q  we  expect  a  tendency  toward  cation 
ordering  above  and  compound  formation  below  Tc. 

10.1  RECENT  ALLOY  FORMATION  MODELS 

There  have  recently  been  several  calculations  of  a  for 
III-V  alloys  based  on  the  idea  that  most  or  all  of  the  heat 
of  mixing  is  due  to  bond  distortion.  Fedders  and  Muller  [1] 
calculate  n  under  the  assumption  that  it  is  due  to  the  changes 
in  the  bond  lengths .  a  and  b,  of  the  binary  compounds,  AC  and 
BC,  required  to  form  a  virtual  crystal  lattice.  This  gives 


n  =  a 


>a-bt2 

a+b 


f 


where  A  can  be  determined  from  the  elastic  properties  of  the 
material.  Mikkelsen  and  Boyce  [2]  have  found  from  EXAFS 
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measurements,  however,  that  the  AC  and  BC  bonds  in  alloys 
are  within  20  percent  of  their  values,  a  and  b,  in  the  binary 
compounds,  which  reduces  X  by  about  20. 

To  maintain  essentially  binary  bond  lengths,  a  and  b, 
in  a  virtual  crystal  alloy  the  bond  angles  must  also  be  distorted 
so  that  the  factor  of  20  reduction  above  is  an  under  estimate 
of  X.  Mikkelson  [3]  has  taken  into  account  both  bond  length 
and  angle  distortions  to  obtain  values  of  n  in  good  agreement 
with  experimental  values  of  0  from  pseudo-binary  phase  diagrams. 
Although  Muller  [4]  has  shown  that  bond  distortion  should 
reduce  disorder  below  that  expected  for  a  random  cation  distribu¬ 
tion  of  A  and  B,  all  elastic  models  necessarily  result  in 
a  positive  heat  of  mixing  since  strain  energy  is  positive 
definite. 

Zunger  et  al.  [5]  have  included  charge  transfer  effects 
with  bond  distortion  in  their  model  of  semiconductor  alloy 
formation.  They  find  that,  although  the  negative  excess 
energy  from  charge  transfer  is  small  in  comparison  to  the 
positive  excess  heat  of  mixing  from  bond  distortion,  the 
chemical  charge  transfer  provides  a  means  of  stabilizing 
long  range  order  and  compound  formation.  This  is  because 
ordered  alloys  or  compounds  can  accomodate  bond  distortion 
with  less  strain  energy  than  disordered  alloys  with  large 
strain  energy  associated  with  random  cation  distribution. 
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As  an  example,  for  the  alloy  InQ  5GaQ  5P  they  obtain 
Q  -  +0 . llev/atom-pair  which  is  the  same  as  obtained  by  Mikkelsen. 
For  the  compound  InGaP2  Zunger  et  al.  obtain  an  equivalent 
excess  energy  parameter  n£  =  -  0.047eV/atom-pair.  On  this  basis 
we  can  qualitatively  compare  alloys  with  In  and  Ga  cations  to 
those  with  Al  and  Ga  cations.  That  is,  alloys  with  In  and  Ga 
cations  are  expected  to  have  relatively  large  bond  distortion 
(atomic  size  differences)  and  large  charge  transfer  (electro¬ 
negativity  differences)  compared  to  alloys  with  Al  and  Ga  cations 
Thus,  the  small  positive  or  zero  heat  of  mixing  in  alloys  with 
Al  and  Ga  cations  is  proportionally  offset  by  small  negative 
or  zero  excess  energy  due  to  chemical  interaction. 

10.2  X-RAY  MEASUREMENTS 

As  mentioned  in  our  last  annual  report  [6]  we  have  obtained 
some  evidence  for  long  range  ordering  in  InxGai_xp  grown  on  GaAs 
by  liquid  phase  epitaxy.  This  is  indicated  in  Figure  1  where 
we  show  back-reflection  Laue  photographs  for  (a)  a  (211) Ga 
oriented  GaAs  substrate  and  (b)  an  InQ  49GaQ  51P  lattice-matched 
epitaxial  layer  using  a  Cu  target.  In  these  measurements 
additional  spots  in  the  photograph  for  the  epitaxial  layer 
could  easily  be  interpreted  as  evidence  for  long  range  order. 

What  we  see  instead  are  enhanced  intensity  0  46  and  0  6  4  spots. 
These  enhanced  intensity  spots  were  also  observed  in  about 
eight  other  InxGa]_.xp  layers  with  x  between  0.49  and  0.52,  but 
were  not  observed  in  one  A1q  4GaQ  ^As  layer  we  measured.  When 
Laue  photographs  f or  an  InQ  4gGaQ  51P  layer  were  taken  with  a 


Figure  1  Back-reflection  Laue  photographs  for  (a)  a  {211} Ga  orientated 
GaAs  substrate  and  (b)  an  Ino . 49Ga0 . 51p  lattice-matched 
epitaxial  layer.  The  primary  X-ray  beam  from  a  Cu  target 
is  in  the  [ Oil J  direction. 
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Mo  target ,  the  0  4  6  and  0  6  4  spots  were  not  enhanced  but 
additional  2,10  12,  2  10,12,,  2  12,10,,  and  2  12  10  were  observed. 


Table  I.  Reflections  which  satisfy  the  Bragg  condition  for 
characteristic  Ka  lines  near  the  lattice  constant 
of  m0,5Ga0,5P. 


Equivalent 

Cubic 

Reflection 

Structure  Factors 

Characteristic 

Bragg 

Condition 

Sphalerite 

Chalcopyrite 

Layered 

032 

0 

0 

fIn"fGa 

2nd  Order 

Cu  Ka 

064 

f  In  +  fGa“2fp 

fIn+fGa”2fp 

fIn+fGa"2fp 

1st  Order 

Cu  Ka 

165 

0 

0 

fIn“fGa 

2nd  Order 

Mo  Ka 

f In  +  fGa*2fp 

£In+£Ga+2£p 

£Xn+£Ga+2£p 

1st  Order 

Mo  Ka 

i - 

2127 

0 

fIn“fGa 

0 

1st  Order 

Cu  Ka 

4109 

0 

f In”fGa 

0 

1st  Order 

Cu  Ka 

In  Table  I  we  compare  the  structure  factors  for  several 
reflections  in  the  sphalerite  structure  and  two  of  the  most 
likely  ordered  structures  for  InGaP2,  chalcopyrite  (space 
group  14  2d)  and  primitive  layered  (space  group  P4  2m)  structures. 
First,  note  that  superlattice  structure  factors  are  small 
(fjn“fGa)  and  will  be  observed  only  when  they  satisfy  the 


V* 


Bragg  condition,  2dn^^  sin  6  =  n  A ,  for  one  of  the  high  intensity 
characteristic  x-ray  lines.  For  example,  the  2  12  7  and  4  10 9 
(equivalent  cubic  notation)  chalcopyrite  reflections  satisfy 
the  Bragg  condition  in  1st  order  for  Cu  Ka  (A  =  1.5418A)  lines 
near  the  lattice  constant  of  GaAs  (5.6532A).  Since  we  easily 
observe  these  reflections  in  our  chalcopyrite  ZnSnP2  (5.6407A) 
epitaxially  grown  on  GaAs,  we  conclude  that  the  In  Ga,  P  is 
not  ordering  in  a  chalcopyrite  structure. 

Second,  notice  that  the  064  set  of  reflections  are  allowed 
in  all  three  crystal  structures  and  satisfy  the  Bragg  condition 
in  1st  order  for  Cu  Ka.  If  the  InxGa1_xP  has  a  layered 
structure,  these  064  reflections  could  be  enhanced  by  032 
reflections,  which  satisfy  the  Bragg  condition  in  2nd  order 
for  Cu  Ka.  In  a  similar  manner  the  2JJ2  10  reflections,  which 
satisfy  the  Bragg  condition  in  1st  order  for  Mo  Ka(A=  0.7107A) , 
could  be  enhanced  by  165  reflections,  which  satisfy  the  Bragg 
condition  in  2nd  order  for  Mo  Ka . 

Using  tabulated  atomic  scattering  factors  [7]  we  can 
calculate  the  intensity  of  the  064  Laue  spots  that  would 
be  expected  with  different  materials  and  crystal  structures. 
Taking  into  account  032,  064,  096,  and  0  12.8  reflections, 
the  following  reflected  intensity  ratios  are  obtained. 

InGaP2  InQ  5GaQ  5P  InGaP2 

InQ  5Ga 0  5P  GaAs  GaAs 

3.5  60  220 

Although  this  calculation  neglects  absorption  and  secondary 

extinction,  it  shows  that  the  064  enhancement  shown  in 


-50- 


Figure  1  could  be  simply  due  to  the  difference  between 
sphalerite  InQ  5GaQ  5P  and  GaAs  and  not  layered  InGaP2 
and  GaAs. 

10.3  LAYERED  STRUCTURE 

Recently,  Kuan,  Kuech  and  co-workers C 8 ]  have  reported 
long-range  order  in  A1  Ga,  As.  They  claim  to  see  the  most 
pronounced  evidence  for  ordering  from  A1q  75Gag  25As  grown 
on  {110} GaAs  at  7008C  by  organo-metallic  vapor  phase  epitaxy. 
Electron  diffraction  measurements,  which  are  sensitive  to 
surface  structure,  were  used  to  indicate  the  presence  of 
a  layered  crystal  structure  (the  P4  2m  space  group  discussed 
above)  with  c-axes  orientations  in  the  growth  plane.  No 
significantly  different  physical  properties,  however,  were 
obtained  for  this  material. 


A  fairly  definitive  x-ray  experiment  which  would  indicate 
the  presence  of  a  layered  crystal  structure  is  shown  schematic¬ 
ally  in  Figure  2.  The  geometry  of  the  experiment  is  the 
same  as  that  used  to  obtain  the  back-reflection  Laue  photographs 
shown  in  Figure  1.  The  layered  face  centered  cubic  cation 
sublattice  with  cations  as  open  and  closed  circles  and  [001] 
c-axes  is  shown  as  observed  in  the  [100]  direction.  Figure  2  (a) 
shows  that  for  this  c-axis  orientation  reflections  from  the 
(046)  planes  between  the  (023)  planes  do  not  cancel  out  (023) 
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Figure 


2  Schematic  of  a  back- ref lection  Laue  measurement 
which  would  indicate  the  presence  of  a  layered 
crystal  structure.  The  orientation  of  the 
experiment  is  the  same  as  that  used  to  obtain 
the  photographs  in  Figure  1. 


reflections.  Thus,  both  046  and  02  3  superimpose  on  the 
bottom  of  the  film.  Figure  2  (b)  shows,  however,  that  for 
the  same  c-axis  orientation  reflections  from  the  (064)  planes 
cancel  out  reflections  from  the  (032)  planes.  Thus,  only 
06  4  is  observed  on  the  top  of  the  film.  A  significant 
difference  in  intensity  of  these  spots  on  the  film  (we  calcu¬ 
late  3.5  for  In  and  Ga  cations),  therefore,  would  indicate 
a  layered  structure  with  one  c-axis. 

If  the  c-axis  in  Figure  2  were  in  the  [010]  direction, 
one  can  see  that  the  0  6  4  and  0  3  2  at  the  top  of  the  film 
would  superimpose  while  the  023  reflection  would  be  cancelled 
by  0  64  at  the  bottom.  For  the  c-axis  in  the  [100]  direction 
both  0  3  2  and  0  23  reflections  would  be  cancelled,  since 
the  layered  structure  looks  like  a  sphalerite  structure  to 
an  [Oil]  x-ray  beam  orthogonal  to  the  c-axis.  The  examination 
of  a  tetragonal  structure  with  fixed  primary  x-ray  beam 
orientation  and  three  c-axis  orientations  is  equivalent  to 
the  examination  of  a  tetragonal  structure  with  fixed  c-axis 
orientation  and  three  primary  x-ray  beam  orientations.  The 
latter,  of  course,  is  experimentally  achievable  and  we 
summarize  potential  observations  in  Table  II. 


Table  II  Summary  of  conclusions  that  can  be  obtained 
from  Laue  measurements  with  primary  beam  in 
three  < 110  >  directions . 


Beam 

Directions 

Enhanced 

0  6  4  Spots 

P4  2m 
c-axes 

[Oil] 

0  4  6 

— 

[101] 

— 

40l 

[001] 

[110] 

— 

— 

[Oil] 

04  6 

06  4 

[001] 

[101] 

— 

406 

and 

[110] 

4  6  0 

— 

[010] 

[Oil] 

046 

06  4 

[001] 

[010] 

[101] 

WMW 

406 

and 

[110] 

640 

460 

[100] 

The  salient  feature  of  the  experiments  indicated  in 
Table  II  is  that,  to  determine  the  number  of  c-axes  in  a 
layered  sample,  x-ray  measurements  must  be  taken  with  the 
primary  beam  in  at  least  two  of  the  <110>  directions 
indicated.  At  the  present  time  we  have  looked  at  our 
InxGa1_xP  samples  with  the  primary  beam  in  only  one  <110> 
direction.  For  these  samples  no  asymmetry  in  064  spot 
intensity  has  been  observed  which  could  not  be  attributed 
to  slight  misorientation  of  the  primary  beam. 
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10.4  CONCLUSIONS 

From  the  discussion  above  we  conclude  that  charge  transfer 
effects  in  III-V  alloys  may  provide  a  mechanism  for  stabilizing 
equilibrium  long  range  order  and  compound  formation.  The 
resulting  total  negative  excess  energy  involved,  however, 
appears  to  be  sufficiently  small  that  low  growth  temperatures 
and  kinetic  effects,  such  as  substrate  orientation,  may  be 
critical  factors.  This  appears  to  be  the  case  for  the  recently 

reported  ordering  in  AlxGa1_xAs  and  should  also  be  true  for  alloys 
with  larger  atomic  size  and  electronegativity  differences.  We  are 
now  in  the  process  of  constructing  an  OMVPE  reactor  to  examine 
these  effects. 
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